Topology of quantum vacuum

Topology of quantum vacuum

G.E. Volovik

Abstract Topology in momentum space is the main characteristics of the ground
states of a system at zero temperature, the quantum vacua. The gaplessness of
fermions in bulk, on the surface or inside the vortex core is protected by topology,
' and is not sensitive to details of the microscopic physics (atomic or trans-Planckian).

Irrespective of the deformation of the parameters of the microscopic theory, the en-
^ I ergy spectrum of these fermions remains strictly gapless. This solves the main hier-

. archy problem in particle physics: for fermionic vacua with Fermi points the masses

of elementary particles are naturally small. The quantum vacuum of Standard Model
is one of the representatives of topological matter alongside with topological super-
fluids and superconductors, topological insulators and semi-metals, etc. There is a
number of of topological invariants in momentum space of different dimensions.
They determine universality classes of the topological matter and the type of the
effective theory which emerges at low energy. In many cases they also give rise to
emergent symmetries, including the effective Lorentz invariance, and emergent phe-
nomena such as effective gauge and gravitational fields. The topological invariants
in extended momentum and coordinate space determine the bulk-surface and bulk-
vortex correspondence. They connect the momentum space topology in bulk with
the real space. These invariants determine the gapless fermions living on the surface
of a system or in the core of topological defects (vortices, strings, domain walls,
solitons, monopoles, etc.). The momentum space topology gives some lessons for
quantum gravity. In effective gravity emerging at low energy, the collective variables
are the tetrad field and spin connections, while the metric is the composite object of
tetrad field. This suggests that the Einstein-Cartan-Sciama-Kibble theory with tor-
sion field is more relevant. There are also several scenarios of Lorentz invariance
violation governed by topology, including splitting of Fermi point and development
of the Dirac points with quadratic and cubic spectrum. The latter leads to the natural
emergence of the Hof ava-Lifshitz gravity.

Low Temperature Laboratory, Aalto University, P.O. Box 15100, FL00076 AALTO, Finland
L.D. Landau Institute for Theoretical Physics, Kosygina 2, 1 19334 Moscow, Russia
e-mail: vo lovik@boo jum . hut ■ £i |

1

2

G.E. Volovik

1 Introduction

Topological approach seeks to recover the Standard Model of particle physics and
general relativity as the low-energy phenomena emerging in the complicated micro-
scopic media called the quantum vacuum. Momentum space topology distributes
the quantum vacua into universality classes, which are characterized by different
types of topological invariants, expressed in terms of propagators - Green's func-
tions. The similar invariants characterize the ground states of topological materials:
topological superfluids, superconductors, insulators, semi-metals, etc. The only dif-
ference between the topological materials and the quantum vacuum of the Standard
Model is that the latter obeys several symmetries, which are absent in condensed
matter physics. These are the Lorentz invariance and different types of gauge sym-
metry. However, in principle, the symmetries of the Standard Model can arise as
the emergent phenomena. This is precisely what happens in the universality class of
vacua, which are characterized by the existence of the topologically protected Fermi
point in momentum space - the point where the energy of the fermionic excitations
is necessarily nullified due to non-zero topological invariant. Such point exists in
the vacuum of the Standard Model in the symmetric phase above the electroweak
interaction, and also in some topological materials, whose ground states also be-
long to the Fermi point universality class. These are superfluid ^He in the phase A,
and topological semimetals. That is why we call the symmetric phase of the Stan-
dard Model as quantum vacuum in the semi-metal state. At low energy, i.e. in the
vicinity of the point node in the spectrum of fermionic excitations, these excitations
behave as relativistic particles - chiral (left handed or right handed) Weyl fermions
- interacting with effective gauge fields and gravitational fields. In principle, all the
ingredients of Standard Model and general relativity can be recovered as phenomena
naturally emerging in the background of the vacuum (ground state), which belongs
to the Fermi point universaUty class.

Here we consider also the other universaUty classes of topological media, illus-
trating them on condensed matter and particle physics examples. Among them is the
vacuum of Standard Model below the electroweak transition, where the momentum-
space topology prescribes the appearance of masses for all fermionic excitations:
quarks and leptons. This massive phase of the Standard Model vacuum is also topol-
gically nontrivial: it belongs to universality class, which includes such topological
media as superfluid ^He in the phase B and three-dimensional topological insula-
tors obeying time reversal symmetry. The fermionic spectrum in all these systems is
fully gapped, but there are topologically protected gapless (massless) states on the
surface of these media or at the interfaces. That is why we call the massive phase of
the Standard Model as the quantum vacuum in the insulating state.

Another interesting class of topological media includes the vacua where the
fermionic excitations are exotic gapless Dirac particles with non-relativistic quadratic,
cubic or higher order spectrum. The effective quantum electrodynamics, which
emerges in such vacua, experiences the so-called anisotropic scaUng, in which the
space and time transform in different ways. The topological mechanism of emergent

Topology of quantum vacuum

3

anisotropic scaling can serve as a source of the recently suggested quantum gravity
at short distances (the so-called Hofava-Lifshitz gravity).

Some sections are devoted to the fermionic bound states leaving on the bound-
ary of the topological media, on the interfaces between the vacua with different
topological invariants, or inside the topological defects, such as strings in relativis-
tic theories and vortices in superfluids and superconductors. The topology can say
may or may not these bound states have exactly zero energy. In superfluids and
supercondictors such zero-energy states behave as Majorana fermions, which are
still elusive in particle physics. The topological consideration of bound states sug-
gests another scenario of emergence of the chiral fermions of the Standard Model:
they may emerge as the fermion zero modes at the interface between two different
insulating quantum vacua in the 4-+1 space-time.

2 Quantum vacuum as topological medium

2.1 Symmetry vs topology

There is a fundamental interplay of symmetry and topology in physics, both in con-
densed matter and relativistic quantum fields. Traditionally the main role was played
by symmetry: gauge symmetry of Standard Model and GUT; symmetry classifica-
tion of condensed matter systems such as solid and liquid crystals, magnets, super-
conductors and superfluids; universality classes of spontaneously broken symmetry
phase transitions; etc. The last decades demonstrated the opposite tendency in which
topology is becoming primary being the main characteristics of quantum vacua in
relativistic quantum fields and their condensed-matter counterparts - ground states
of the condensed-matter systems at T = 0, see reviews |[Tl|2l[3]|4l and earlier papers

13 |6l El 111 IE [13 III] El [!3 [II [B] [111 [nl [11] .

Topology describes the properties of a system, which are insensitive to the details
of the microscopic physics (atomic physics in condensed matter and Planck-scale
physics of quantum vacuum). The momentum-space topology, which we discuss in
this Chapter, determines universality classes of the topological media and the ef-
fective quantum field theory which emerges in such media at low energy and low
temperature, including the type of the energy spectrum of fermionic excitations (in
condensed matter systems fermionic excitations play the role of elementary parti-
cles). The topology also gives rise to emergent symmetry, the symmetry which did
not exist on a microscopic level. Examples are provided by the point nodes (zeroes)
in the energy spectrum of fermionic excitations: some of the nodes are protected
by topology, i.e. they are robust to the deformations of the system. Close to such
nodes the spectrum of fermions becomes "relativistic", i.e. the spectrum forms the
Dirac cone, and fermionic excitations behave as Weyl, Dirac or Majorana particles.
As a result the effective relativistic quantum field theory emerges, in which bosonic
collective modes give rise to effective gauge field and effective metric. All this is

4

G.E. Volovik

the consequence of the topological theorem - the Atiyah-Bott-Shapiro construction

im.

Among the existing and potential representatives of topological materials one can
find those in which the spectrum of fermionic excitations has a gap (the gap corre-
sponds to mass of a relativistic particle). These systems include 3D topological band
insulators |3l; fully gapped superluid ^He-B lfT4l [T9l l20l : 2D materials exhibiting
intrinsic (i.e. without external magnetic field) quantum Hall and spin-Hall effects,
such as gapped graphene ifTsl : thin film of superluid ^He-A and quasi 2D planar
phase of triplet superfluid lfT6l 1211 [iTl ; and chiral superconductor Sr2Ru04 ll22ll .
These materials have the topological properties similar to that of the quantum vac-
uum of Standard Model in its massive phase, i.e. in the state below the electroweak
transition, where all elementary particles are massive |[23l. The latter state can be
called the "insulating" state of the relativistic quantum vacuum.

The gapless topological media are represented by superluid ^He-A HI; topo-
logical semimetals EH |25] EH |27| ; gapless graphene 118] [lH 12 EO] EH ; nodal
cuprate superconductors ll28l : and noncentrosymmetric |[32l, l33l superconductors.
These materials are similar to the quantum vacuum of Standard Model in the state
above the electroweak transition, where all elementary particles are massless lfn i23l .
This state can be called the "semimetal" phase of the relativistic quantum vacuum.

2.2 Green's function vs order parameter

Topology operates in particular with integer numbers - topological invariants or
topological charges - which do not change under small deformation of the system.
The topological invariants appear in very different situations. The commonly known
are topological charges, which describe the topological defects - the inhomogeneous
objects such as quantized vortices, dislocations, domain walls, hedgehogs, solitons,
etc. in condensed matter systems or the cosmic strings, magnetic monopoles, instan-
tons, etc. in particle physics. The topological charge protects the topological defects
from destruction, while the summation law of the topological charges regulates pro-
cesses of scattering, merging, splitting, and other transformations of these defects in
their dynamics.

Here we consider the topological invariants of different type: they describe the
homogeneous ground states or the homogeneous quantum vacua. The conservation
of these topological charges in particular protects from destruction different types of
nodes in the fermionic energy spectrum of a system. The node in the spectrum is an
object in momentum space, and thus we call these invariants as momentum-space
invariants, as distinct from the real-space invariants describing topological defects.
In other words, the real-space invariants describes the topologically nontrivial inho-
mogeneous configurations of the order parameter fields f (r,f) in space-time, while
the momentum-space invariants describe the nontrivial momentum-space config-
uration of the propagator - Green's function G(p, o) or other response function,
which characterizes the homogeneous ground state of a system (the vacuum state)

Topology of quantum vacuum

Fermi surface: vortex ring in p-spa ce

fully gapped topological matter:

metals,
normal 3 He

skymiion in p-space

-0

3He-B, topological insulators,
quantum-Hall states, 3He-A film,
vacuum of Standard Model

Weyl point - hedgehog in p-space
3He-A, vacuum of SM. topological semimetals (Abrikosov)

CO

P 'P

— ---^ ^ "

P

=P2

Fermi arc on

Py

3He-A surface

_ r

Pz f-I^l''

••Tif* Fermi surface ••,'•*'*

t

*

Dirac strings in p-space terminating on monopole

4

Pz .Jt..

flat band ''j" %"*

t

on 3He-A vortex *

Khodel-Shaginyan flat band: Tt-vortex in p-space

Fig. 1 Topological matter, represented in terins of topological objects in moinentum space, (top
lefty. Fermi surface is the momentum-space analog of the vortex line: the phase of the Green's
function changes by Itc around the element of the line in (co,p)-space. (top middle): Fermi point
(Weyl point) is the counterpart of a hedgehog and a magnetic monopole. The hedgehog in this
figure has integer topological charge N = +1, and close to this Fermi point the fermionic quasi-
particles behave as Weyl fermions. Nontrivial topological charges in terms of Green's functions
support the stability of the Fenni surfaces and Weyl points with respect to perturbations including
interactions (11 1341 . In terms of the Berry phase 1361 the Fermi point represents the p-space coun-
teipart of Dirac magnetic monopole with unobservable Dirac string (see Ref . 1371 and Fig. 1 1 .4 in
fTl). (top right): Topological insulators and fully gapped topological superfluids/superconductors
are textures in momentum space: they have no singularities in the Green's function and thus no
nodes in the energy spectrum in the bulk. This figure shows a skyrmion in the two-dimensional
momentum space, which characterizes two-dimensional topological insulators exhibiting intrinsic
quantum Hall or spin-Hall effect, (bottom left): Flat band emerging in strongly interacting systems
1411 . This dispersionless Fermi band is analogous to a soliton terminated by half-quantum vortices:
the phase of the Green's function changes by n around the edge of the flat band 1421 . (bottom right):
Fermi arc on the surface of ^^He-A 1431 and of topological semi-metals with Weyl points |27||26l
and flat band inside the vortex core of^He-Al44l serve as the momentum-space analog of a Dirac
string terminating on a monopole. The Fenni surface formed by the surface bound states terminates
on the points where the spectrum of zero energy states merge with the continuous spectrum in the
bulk, i.e. with the Weyl points.

6

G.E. Volovik

im [T8l |34] [35l . The Green's function is generally a matrix with spin indices. In ad-
dition, it may have the band indices (in the case of electrons in the periodic potential
of crystals).

Due to topological stability, the nodes in the spectrum survive when the interac-
tion between the fermions is introduced and/or modified. Thus the momentum-space
topology is the main reason why there are gapless quasiparticles in condensed mat-
ter and (nearly) massless elementary particles in our Universe. Topology explains
the masslessness of the quantum vacuum in its semimetal phase, and allows to solve
the hierarchy problem in particle physics: if the masslessness (or gaplessness) of the
quantum vacuum is not protected, the natural values of the masses of elementary
particles are on the order of the Planck energy scale, which is in huge disagreement
with masses of elementary particles in the insulating vacuum. However, if masses
appear as a result of the spontaneously broken symmetry in the initially gapless
vacuum, they can be small.

The momentum-space topological invariants are in many respects similar to their
real-space counterparts, which describe topological defects in condensed matter sys-
tems and quantum vacuum (see Fig.[T]i. In particular, the point node in the spectrum
(the Fermi point or Weyl point, see Sec. [3]i is the counterpart of the real-space point
defects, such as hedgehog in ferromagnets or magnetic monopole in particle physics
(Fig.[T]fo/7 center). The summation law of the momentum-space topological charges
regulates processes of scattering, merging, splitting, and other transformations of the
nodes in spectrum during deformation of the vacuum state. Another example is the
Fermi surface in metals, which we consider in Sec. 12.31 It is topologically stable,
because it is analogous to the vortex loop in superfluids or superconductors (Fig.[T]
top left).

The fully gapped topological matter, such as topological insulators and fully
gapped topological superfluids have no nodes in the spectrum and thus no singu-
larities in the Green's function. They are the counterpart of the real-space objects
which do not have singularities, but still are topologically non-trivial - the so-called
textures or skyrmions. Thus the topological insulator and the vacuum of Standard
Model in insulating phase can be considered as the momentum-space skyrmion
(Fig. 1 top right).

The further extension of topology to the combined real + momentum space al-
lows us to consider topologically protected spectrum of the so called fermion zero
modes - fermions living on the real-space topological objects such as domain walls,
strings and monopoles. The fermion zero modes are protected by topological in-
variants expressed in terms of the Green' s function in the extended phase space,
G(co,p;f,r) ||38llIl|39l[I]|40l|3l.

2.3 Fermi surface as topological object

We shall start with gapless vacua - the vacua whose fermionic excitations has no
gap. In the gapless media the fermionic degrees of freedom are not frozen out at

Topology of quantum vacuum

7

CO

p.

C

Fig. 2 Nodes of co-dimension 1 in 2+1 systems. Green's function has singularities on line £0 = 0,
P? + p^, = in the three-dimensional space ((B,p, , p, ). Stability of Fermi surface is protected by
the invariant {TJ which is represented by integral over an arbitrary contour C around the Green's
function singularity. This is applicable to nodes of co-dimension 1 in any D+\ dimension. For
Z) = 3 the nodes form conventional Fenni surface in metals and in normal ^^He.

low temperature. This is very important for practical purposes in electronic devices
and for the existence of life in our Universe. It is the topology which provides the
protection for the gaplessness in condensed matter systems and masslessness of
elementary particles in Standard Model.

For the topological classification of the gapless vacua, the Green's function is
considered on imaginary frequency axis. On the real axis singularities are always
present, for example in the simplest case of non-interacting fermions the poles in
the Green's function, G = 1/ (o — e(p)), are present even in the fully gapped sys-
tem. To distinguish between the gapped and gapless states of the vacuum, and also
between different types of gapless vacua, we should move to the imaginary axis.
Then the simplest Green's function becomes G = 1/ {i(0 — e(p)), which may have
singularity at O) = and e(p) = 0. The latter is possible only for the gapless vacua,
where the energy of excitations e(p) has zeroes. Thus zeroes in the spectrum of
fermionic excitations are manifested as singularities in the Green's function, and we
are interested in such singularities in the Green's function, which are protected by
topology, i.e. are described by topological charge and thus are robust to perturba-
tions of the vacuum state. The transition to imaginary axis thus allows us to consider
only the relevant singularities in the Green's function and to avoid the trivial singu-
larities on the mass shell.

Zeroes in the energy spectrum of fermionic excitations may form different p-
dimensional manifolds in momentum space: there can be zero-dimensional points
p — Q, one-dimensional lines of zeroes p = I, two-dimensional surfaces p ~2 and
even the whole three-dimensional bands with zero energy p ^3. Since the dimen-
sion D of the momentum space can be also different, for topological classification it

8

G.E. Volovik

is more instructive to use the co-dimension of zeroes. By co-dimension we denote
D — p, i.e. the dimension D of p-space minus dimension p of the nodes in energy
spectrum.

We start with zeroes of co-dimension 1 . This refers to two-dimensional Fermi
surface in three-dimensional metal (3 — 2 = 1), D=l Fermi line in D=2 systems
(2—1 = 1) and Fermi point in D = 1 systems (1 — = 1). The general analysis ifTSl
demonstrates that topologically stable nodes of co-dimension 1 are described by the
group Z of integers. The corresponding winding number is expressed analytically
in terms of the Green's function |[T]:

N = tr <f^^.Gico,p)d,G-\co,p) . (1)

Here the integral is taken over an arbitrary contour C around the Green's function
singularity in the D+l momentum-frequency space. See Fig.|2]for D = 2. Example
of the Green's function in any dimension D is scalar function G^'(a),p) = i(0 —
vf(|p| — Pf)- For D = 2, the singularity with winding number = 1 is on the line
(0 = 0, p^ + p^, = p^p, which represents the one-dimensional Fermi surface.

Due to nontrivial topological invariant, Fermi surface survives the perturbative
interaction and exists even in systems without poles in the Green's function and
quasiparticles are not well defined. The systems without poles include the so-called
marginal and Luttinger liquids in condensed matter and the so-called unparticles in
relativistic theories (see Sec. |3.2| l.

3 Vacuum in a semi-metal state

Here we consider the topology of the quantum vacuum of Standard Model in its
symmetric phase above the electroweak transition, where all fermions are massless
Weyl fermions. Standard Model vacuum obeys the relativistic invariance and thus
cannot belong to the universality class of vacua with Fermi surface: the latter vio-
lates the Lorentz symmetry. The gapless (massless) phase of Standard Model can
be consistent with Lorentz symmetry only if the vacuum of Standard Model be-
longs to the universality class of vacua, which is characterized by point nodes in
the fermionic spectrum - the Fermi points. The Fermi point has the co-dimension
3. The Fermi point can be of the chiral type, when the fermionic excitations near
the Fermi point behave as relativistic left-handed or right-handed Weyl fermions,
such as left neutrino, left electron, left quarks and their right-handed counterparts
P31 [l1. In this case the Fermi point is called the Weyl point. The Fermi point can
be also of the Dirac type, when the fermionic excitations behave as massless Dirac
fermions. There is also the class of quantum vacua in which the nodal point obeys
the Z2 topology, with the topological charge obeying the summation rule 1 + 1=0.
In such vacuum the fermionic excitations behave as relativistic massless Majorana
particles ifTsl .

Topology of quantum vacuum

3.1 Fermi points in 3+1 vacua

9

Fermi point is the Green's function singularity described by the following topolog-
ical invariant expressed via integer valued integral over the surface a around the
singular point in the 4-momentum space = (fi),p) HI:

^ = ^ I ^^Pp G ' ^^P^ G ' ^^/^v G '. (2)

If the invaraint (O is nonzero, the Green's function has a singularity inside the sur-
face C7, and this means that fermions are gapless. The typical singularities have
topological charge A^=:+l orA^=— l.In the vicinity of the Fermi point with such
topological invariant the spectrum of fermionic excitations reproduces the spectrum
the chiral Weyl fermions, right-handed or left-handed respectively. This is the con-
sequence of the so-called Atiyah-Bott-Shapiro construction ifTSl . which leads to the
following general form of expansion of the inverse fermionic propagator near the
Fermi point with A^ = +lorA'=— 1:

G-\p,)=e^ar"{pp-pf^) + ---. (3)

Here = (1 , CJj, CJy, o;) are Pauli matrices (or Dirac matrices in the more general
case); the expansion parameters are the vector p^^ ^ indicating the position of the
Fermi point in momentum space where the Green's function has a singularity, and
the matrix ; ellipsis denote higher order terms in expansion. The expansion results
in the linear spectrum of fermionic excitations, which induces the effective Lorentz
invariance.

3.2 Emergent relativistic fermionic matter

By continuously adiabatic deformation one may transform equation (O to the equa-
tion which describes the relativistic Weyl fermions

G-\p^) = i(0+Na-p+--- , N = ±l. (4)

Here the position of the Fermi point is shifted to p^^^ = 0; the matrix is deformed
to unit matrix; and ellipsis denote higher order terms in (O and p. This means that
close to the Fermi point with = + 1 , the low energy fermionic excitations behave
as right handed relativistic particles, while the Fermi point with = — 1 gives rise
to the left handed particles. They obey the following effective Weyl Hamiltonian

(5)

10

G.E. Volovik

Fig. 3 Weyl point as the conical point of level crossing. The negative energy branch occupied by
the fermionic excitations touches the positive energy branch.

The spectrum of this Hamihonian is shown in Fig. [3], which demonstrates that
the Fermi or Weyl point is an example of the exceptional point of level crossing
analyzed by von Neumann and Wigner P6l . This analysis demonstrates that two
branches of spectrum, which have the same symmetry, may touch each other at the
conical (or diabolical) point in the three-dimensional space of parameters. In other
words the twofold degeneracy has co-dimension = 3. In our case the two branches
correspond to positive and negative energy states, while the three parameters are the
components of momentum p,, py and |[8]|32l. Topological invariants for points
at which the branches of spectrum merge were introduced by Novikov 1471 .

However, this is not the whole story. Topological invariant of the Weyl point does
not change if the Green's function is multiplied by any power of p^ + o^, which
violates the simple pole structure of the propagator in Eq.(|4li. Such violation may
occur due to the infrared divergences, which take place in quantum field theories. In
this case in the vicinity of Fermi point one has

with 7 7^ 1 . In particle physics, the hypothetical particles described by the modified
Green's function with 77^ 1 are called unparticles 1481 |49l . This modification does
not change the topology of the propagator: the topological charge of singularity is
for arbitrary parameter 7 lO . For fermionic unparticles one has 7 = 5 /2 — t/y, where
du is the scale dimension of the quantum field.

The main property of the vacua with Weyl points is that according to (|4|i, close
to the Weyl points the massless relativistic fermions emerge. This is consistent with
the fermionic content of our Universe, where all the elementary particles - left-
handed and right-handed quarks and leptons - are Weyl fermions. Such a coinci-
dence demonstrates that the vacuum of Standard Model is the topological medium

Topology of quantum vacuum

11

of special universality - the class of vacua which have topologically protected Fermi
points. This solves the hierarchy problem, since the value of the masses of elemen-
tary particles in the vacua of this universality class is identically zero.

Let us suppose for a moment, that there is no topological invariant which protects
massless fermions, and the vacuum of our Universe belongs to the class of the fully
gappedvacua. Then the natural masses of fermions must be on the order of Planck
energy scale: M ^Ep ^ lO'^ GeV. In such a natural Universe, where all masses are
of order E-p, all fermionic degrees of freedom are completely frozen out because of
the Bolzmann factor e^'^l'^ , which is about at the temperature corresponding

to the highest energy reached in accelerators. There is no fermionic matter in such
a Universe at low energy. That we survive in our Universe is not the result of the
anthropic principle (the latter chooses the Universes which are fine-tuned for life
but have an extremely low probability). Our Universe is also natural and its vacuum
is generic, but it belongs to the universality class of vacua with Fermi points. In
such vacua the masslessness of fermions is protected by topology (combined with
symmetry, see below).

3.3 Emergent gauge fields

The vacua with Fermi-point suggest a particular mechanism for emergent symme-
try. The Lorentz symmetry is simply the result of the linear expansion: this sym-
metry becomes better and better when the Fermi point is approached and the non-
relativistic higher order terms in Eq.© may be neglected. This expansion demon-
strates the emergence of the relativistic spin, which is described by the Pauli matri-
ces. It also demonstrates how gauge fields and gravity emerge together with chiral
fermions. The expansion parameters z?^"' and may depend on the space and time
coordinates and they actually represent collective dynamic bosonic fields in the vac-
uum with Fermi point. The vector field p^*'^ in the expansion plays the role of the
effective U{\) gauge field acting on fermions.

For the more complicated Fermi points with |A^| > 1 the shift p^'^ becomes the
matrix field; it gives rise to effective non-Abelian (Yang-Mills) SU {N) gauge fields
emerging in the vicinity of Fermi point, i.e. at low energy IH. For example, the
Fermi point with N = 2 may give rise to the effective SU{2) gauge field in addition
to the effective U{1) gauge field

G'\Pn)=ear" {pji -giAji -g2A|3 -t) + higher order terms , (7)
where t are Pauli matrices corresponding to the emergent isotopic spin.

12

G.E. Volovik

3.4 Emergent gravity

The matrix field in (|7|i acts on the (quasi)particles as the field of vierbein, and thus
describes the emergent dynamical gravity field. As a result, close to the Fermi point,
matter fields (all ingredients of Standard Model: chiral fermions and Abelian and
non-Abelian gauge fields) emerge together with geometry, relativistic spin, Dirac
matrices, and physical laws: Lorentz and gauge invariance, equivalence principle,
etc. This is the result of the natural coarse graining: in the low energy corner, when
the Fermi point is approached, the huge number of microscopic (Planckian) degrees
of freedom of the physical quantum vacuum is reduced to rather few low-energy
degrees of freedom (chiral fermions, gauge fields and gravity field). All other de-
grees of freedom are highly massive and are completely frozen out at low energy
and temperature.

In such vacua, gravity emerges together with matter If this Fermi point mech-
anism of emergence of physical laws works for our Universe, then the so-called
"quantum gravity" does not exist. The gravitational degrees of freedom become
separated from all other degrees of freedom of quantum vacuum only at low energy.

In this scenario, classical gravity is a natural macroscopic phenomenon emerging
in the low-energy corner of the microscopic quantum vacuum, i.e. it is a typical and
actually inevitable consequence of the coarse graining procedure discussed above.
It is possible to quantize gravitational waves to obtain their quanta - gravitons, since
in the low energy corner the results of microscopic and effective theories coincide.
It is also possible to obtain some (but not all) quantum corrections to Einstein equa-
tion and to extend classical gravity to the semiclassical level. But one cannot obtain
"quantum gravity" (i.e. microscopic physics of the vacuum) by quantization of Ein-
stein equations for the gravitational degrees of freedom, since all other degrees of
freedom of the quantum vacuum will be missed in this procedure.

3.5 Topological invariant protected by symmetry in Standard Model

Standard Model contains equal number hr ~ ni of right and left Weyl fermions, 8
right and 8 left in each generation: «r = «l = 8«g, where iig is the number of genera-
tions (we do not consider Standard Model with Majorana fermions, and assume that
in the insulating state of Standard Model neutrinos are Dirac fermions). For such
Standard Model the total topological charge, obtained after trace is taken over all
fermionic flavors in invariant (|2]i, vanishes, = 8;?;, — 8/1^, = 0. Thus the invariant
^ cannot protect nodes in the spectrum of Standard Model fermions. This would
mean that any interaction between the fermions may make them massive.

However, in the symmetric phase of Standard Modle above the electroweak tran-
sition this does not happen. This is because there is another topological invariant,
which takes into account the symmetry of the vacuum. The gapless state of the vac-
uum with N ^ can be protected by the following topological invariant which is

Topology of quantum vacuum

13

supported by symmetry |[T]:

(8)

K I dS" Gdp^ G ' Gdp^ G- ' Gd,,^ G-

where Kjj is the matrix of symmetry transformation, which either commutes or anti-
commutes with the Green's function matrix. In the symmetric phase of Standard
Model there are two relevant symmetries, both are the Z2 groups, ~ \. One
of them is the center subgroup of SU {2)l gauge group of weak rotations of left
fermions, where the element K is the gauge rotation by angle 2k, K ~ e'^'^^'-. The
other one is the group of the hypercharge rotation be angle 6n, K = e'^^^ . In the
G(224) Pati-Salam extension of the G(213) group of Standard Model, this symme-
try comes as combination of the Z2 center group of the SU {2)r gauge group for
right fermions, e"^^3«, and the element e^Tr^B-i) of jjjg ^4 center group of the SU{4)
color group - the Pm parity (on the importance of the discrete groups in particle
physics see ifSOl |5T1 and references therein). Each of these two Z2 symmetry op-
erations changes sign of left spinor, but does not influence the right particles. Thus
these matrices are diagonal, Kjj = diag(l, 1,..., — 1, — 1,...), with eigen values 1 for
right fermions and —1 for left fermions.

In the symmetric phase of Standard Model, both matrices commute with the
Green's function matrix Gij, as a result Nk in (O is topological invariant: it is robust
to deformations of Green's function which preserve the symmetry K. The value of
this invariant A^/f = 16«g, which means that all 16«g fermions are massless.

Topological invariant protected by symmetry is responsible also for the effect of
chiral anomaly - the anomalous nonconservation of a chiral current ||52ll53l . This
anomaly can be the source of the production of baryonic charge B in the early Uni-
verse, which leads to the excess of matter over anti-matter in the present Universe
(the so-called electroweak baryogenesis, see the recent review ll54l ). The baryo-
production can be expressed in terms of the invariant Nk- in case of the creation of
the baryonic charge by the hypermagnetic U{1) field the relevant symmetry opera-
tor K in Eq.® is A" = BY^, where Y is the generator of symmetry group. The
creation rate of baryons in the hypermagnetic U{1) field is

B=^By-Ey , K^BY\ (9)

In condensed matter, the chiral anomaly generated by Weyl fermions in a crystal has
been considered in Ref. Q, while the analog of the baryoproduction in Eq. (|9]l has
been experimentally tested in the superfluid with Weyl points - the ^He-A |

14

G.E. Volovik

marginal Fermi point

iV = -1+1 =

Fermi points witli opposite N marginal Fermi point splits
annihilate each other & form into two separate

massive Dirac fermions topologically protected Fermi points

Fig. 4 (top) In Standard Model the Fermi points with positive N = +i and negative N = — I
topological charges are at the same point p = 0, fonning the marginal Fermi point with N = 0.
Symmetry K between the Fermi points prevents their mutual annihilation giving rise to the topo-
logical invariant (8) with Nk = 2. (bottom left): If symmetry K is violated or spontaneously broken,
Fermi points annihilate each other and Dirac mass is formed, {bottom right): If Lorentz invariance
is violated or spontaneously broken, the marginal Fermi point splits |56| . The topological quantum
phase transition between the state with Dirac mass and the state with splitted Dirac points have
been observed in cold Fermi gas 1571 .

3.6 Higgs mechanism vs splitting of Fermi points

The gapless vacuum of Standard Model is supported by combined action of topol-
ogy and symmetry K, and also by the Lorentz invariance which keeps all the Fermi
points at p = 0.

Explicit violation or spontaneous breaking of one of these symmetries transforms
the vacuum of the Standard Model into one of the two possible vacua. If, for exam-
ple, the K symmetry is broken, the invariant (O supported by this symmetry ceases
to exist, and the Fermi point disappears. All \6n^ fermions become massive (Fig.|4]
bottom left). This is assumed to happen below the symmetry breaking electroweak
transition caused by Higgs mechanism where quarks and charged leptons acquire
the Dirac masses.

If, on the other hand, the Lorentz symmetry is violated, the marginal Fermi point
splits into topologically stable Fermi points with non-zero invariant A^, which pro-

Topology of quantum vacuum 15

Fig. 5 Sketch of topologically protected point nodes in momentum space with topological charge
A' = ±1 in Eq.([2ll in superconfuctors/superfluids of 0{D2) symmetry class |59i . Chiral fermions
emerge in the vicinity of each Fermi points. They have identical emergent Minkowski metric,
hut different orientations of dreibein. The simplest realization of dreibein for each of eight chiral
fermions is shown by aiTows. The vierbein orientations in Ol^Di) symmetry class superconductors
are more complicated than in Figure: one of the vectors in each vierbein is along the cube main di-
agonal, see Eq.in). All four pairs of left and right Weyl fermions have the same quantum numbers,
while their triads can be transformed to each other by rotations and reflection. They are analogous
to generations of Standard Model fermions.

tects massless chiral fermions {Vig.^bottom right). Since the invariant does not
depend on symmetry, the further symmetry breaking cannot destroy the nodes. One
can speculate that in the Standard Model the latter may happen with the electri-
cally neutral leptons, the neutrinos ||56ll . Most interestingly, Fermi-point splitting of
neutrinos may provide a new source of T and CP violation in the leptonic sector,
which may be relevant for the creation of the observed cosmic matter-antimatter
asymmetry ll58l .

3.7 Splitting of Fermi points and problem of generations

An example of the multiple splitting is provided by the model Hamiltonian for
fermions in superconfuctors/superfluids in the state which belongs to 0{D2) sym-
metry class 1591 . The Bogoliubov - de Gennes (BdG) Hamiltonian for fermions in
this spin-singlet tZ-wave superconductor has the form

H = i=(/ _ pI)^, + ii2pl-pl-p])x, + -;,2)T2 . (10)

At < the energy spectrum is fully gapped, for = the node in the spectrum
appears at p = which at > splits into 8 Fermi points at the vertices of cube in

16

G.E. Volovik

momentum space (see Fig. 13.71 ):

p(")^Z£(±x±y±z) , „=!,... ,8. (11)

These nodes have topological charges = ±1 in Eq.(|2]i, and as a result, close to
each of 8 nodes the Hamiltonian is reduced to the Hamiltonian describing Weyl
fermions, 4 left and 4 right:

//W =e<"' • (p-pW)Ti +e("' • (p-pW)T2 + e("^ . (p-pW)T3. (12)

Each Weyl fermion has its own triad (dreibein). Choosing for simphcity pp ~ y/3
one has

= V2(±x±y±z), (13)

e<"^-±2xTyTz, (14)

e<"' = V3(±yTz). (15)

All triads can be transformed to each other by rotations and/or reflection. So in this
model one obtains four identical copies of right and left relativistic Weyl fermions.
They may be considered as analogs of generations of Standard Model fermions,
but with rig = 4. Different, but related mechanism for the origin of generations is
suggested in 1601 .

4 Exotic fermions

In many systems (including condensed matter and relativistic quantum vacua), the
Fermi points with elementary charges = ±1 may merge together forming either
the neutral point with = or point with multiple (i.e. |A^| > 1 ||6T| ). In this
case topology and symmetry become equally important, because it is the symme-
try which may stabilize the degenerate node. Example is provided by the Standard
Model of particle physics, where 16 fermions of one generation have degenerate
Dirac point at p = with the trivial total topological charge A^ = 8 — 8 = 0. In the
symmetric phase of Standard Model the nodes in the spectrum survive due to a
discrete symmetry between the fermions and they disappear in the non-symmetric
phase forming the fully gapped vacuum ITI. In case of degenerate Fermi point with
|A^| > 1, situation is more diverse. Depending on symmetry, the interaction between
fermionic flavors may lead to splitting of the multiple Fermi point to elementary
Dirac points ll56l : or gives rise to the essentially non-relativistic energy spectrum
E± (/7 — > 0) — > ±p^, which corresponds to different scaling for space and time in the
infrared: r — > ^r, f — > b^t. The particular case of anisotropic scaling with A^ = 3 was
suggested by Hofava for quantum gravity at short distances, the so-called Hoi'ava-
Lifshitz gravity ll62l 16311641 . while the anisotropic scaling in the infrared was sug-

Topology of quantum vacuum

17

gested in Ref. ||65l . The topology of the multiple Dirac point provides another pos-
sible realization of anisotropic gravity, which is different from the scenario based
on Lifshitz point in the theory of phase transitions ||66l .

4.1 Dirac fermions with quadratic spectrum

The non-linear spectrum arising near the Fermi point with N = 2 has been dis-
cussed for different systems including graphene, double cuprate layer in high-
Tc superconductors, surface states of topological insulators and neutrino physics
II§7l[Tl|2l|29l|68l|69l|70l|7l]|71. The spectrum of (quasi)particles in the vicin-
ity of the doubly degenerate node depends on symmetry. Let us consider zeroes of
co-dimension 2 - point nodes in 2+1 systems. The node with topological charge
= +2 takes place in bilayered graphene. According to general classification ifTSll .
the topology alone cannot protect zeroes of co-dimension 2: the topological invari-
ant takes place only in the presence of a symmetry. In particular, if we restrict con-
sideration only to real (Majorana) fermions, the nodes obey Z2 topology with sum-
mation law 1 + 1=0 ifTSl . To make the multiple Fermi point possible we need an
additional symmetry K which extends the group Z2 to the full group of integers Z.
The relevant symmetry protected topological invariant is 1731 |2l |29l l74l :

N = tr i dl KG(o)=0,p)diG-\co = Q,p)=tr i dlKJ^-\p)di.y/f{p),
Am JC JC

(16)

where C is contour around the Dirac point in 2D momentum space {p.x,Py)', K is
the relevant symmetry operator; G is the Green's function matrix at zero frequency,
which can be used as the effective Hamiltonian, ^(p) = G^' (o = 0,p); the oper-
ator K commutes or anticommutes with the effective Hamiltonian.

Provided the symmetry K is preserved and thus the summation law for takes
place, one finds several scenarios of the behavior of the system with the total topo-
logical charge = +2.

(i) One may have two fermions with the linear Dirac spectrum, with the nodes
being at the same point of momentum space. This occurs if there is some special
symmetry, such as the fundamental Lorentz invariance.

(ii) Exotic massless fermions emerge. In the 2D systems, these are gapless Dirac
fermions with parabolic energy spectrum, which emerge at low energy:

£±(p)«V- (17)
They are described by the following effective Hamiltonian

^bx,A.)= ((^^_°^^)2 + =(p.?-p?)cT,-2p,p,(T2, (18)

18 G.E. Volovik

(a) (b)

Fig. 6 Illustration of summation rule for momentum-space topological invariant. Splitting ofN = 2
point with quadratic dispersion (a) into four Dirac points: A' = 2= l + 1 + 1— 1 (b).

where (7i and (72 are Pauli matrices. The topological charge of the node at the
point px = Py = is given by the Eq.(fT6]l. where the symmetry operator K is rep-
resented by the Pauli matrix 03. With effective Hamiltonian ( fTSI ) one obtains that
the node with the quadratic spectrum has the charge = 2. In 3D systems, the
corresponding fermions with the topological charge = 2 in Eq. (|2]i are the semi-
Dirac fermions, with linear dispersion in one direction and quadratic dispersion in
the other Ig?) :

£±(p)«±ycV?+pl- (19)

(iii) The Weyl point with = +2 may split either into two Weyl points each with
A^ = +1 (see 1561 for the relativistic 3+1 system) or into four Weyl points (three with
N = +1 and one with N = — I, see Fig.|6]l. The effective Hamiltonian for the latter
case is l|75ll76ll :

The energy spectrum of this Hamiltonian has four Dirac point: the node at p =
has the topological charge N = — I, while three nodes at px + ipy = —se^'^'^'l^ with
integer k have charges A^ = + 1 each, so that the summation rule A^=:l + 1 + 1 — 1=2
does hold.

Note that in options (ii) and (iii) the (effective) Lorentz invariance of the Dirac
point is violated. This suggests that the topological mechanism of splitting of the
Dirac point ||56l or of the formation of the nonlinear dispersion ||67l may lead to
the spontaneous breaking of Lorentz invariance in the relativistic quantum vacuum,
which in principle may occur in the neutrino sector of the quantum vacuum 1,77. ,78l

ED.

Let us consider the Fermi point with higher degeneracy, described by the sym-
metry protected topological invariant N > 2.

Topology of quantum vacuum

19

4.2 Dime fermions with cubic spectrum

Let us consider first the case with N = 3 lISOl . Examples are three families of right-
handed Weyl 2-component fermions in particle physics; three cuprate layers in high-
Tc superconductors; three graphene layers, etc. If the Fermi point is topologically
protected, i.e. there is a conserved topological invariant A^, the node in the spectrum
cannot disappear even in the presence of interaction, but it can split into nodes
with elementary charge N ~ I. The splitting can be prevented if there is a symmetry
in play, such as rotational symmetry. Here we provide an example of such symmetry.

-2 -1 1 2

9

Fig. 7 Splitting of Fermi point with N = 3 into three Fermi points with W = 1 in the model
discussed in iSOj .

We consider 3 species (families or flavors) of fermions, each of them being de-
scribed by the invariant = +1 in Eq.(fT6]l and an effective relativistic Hamiltonian
emerging in the vicinity of the Fermi point

(P) = C7 • P = a,p_, + GyPy . (21)

The matrix K = a, anticommutes with the Hamiltonian. This supports the topolog-
ically protected node in spectrum, which is robust to interactions. The position of
the node here is chosen at p = 0:

(22)

20

G.E. Volovik

The total topological charge of three nodes at p = of three fermionic species is
= +3. Let us now introduce matrix elements which mix the fermions. If these
elements violate symmetry K, then the topological invariant cannot be constructed
and point node will be destroyed, so let the matrix elements obey the symmetry K.
For the general case of the matrix elements, but still obeying the symmetry K, the
multiple node will split into 3 or more elementary nodes, obeying the summation
rule: 3 = 1 + 1 + 1 = 1 + 1 + 1 + 1 — 1 = - (see Fig.|7]l. However, in the presence of
some extra symmetry, which prevents splitting, the branch with the cubic spectrum
emerges.

E/g

Fig. 8 Spectrum of the Hamiltonian <23b showing cubic dispersion for the lowest two eigenvalues
around the point p = 0. The spectrum has been calculated with equal coupling strengths gi2 =
Si3 ~ S2i = g, but cubic spectrum characterized by topological charge N = 3 preserves for any
nonzero values of coupling strengths. Spectrum is shown as function of p, at p,. = 0.

Example is provided by the following Hamiltonian ISOl

^(P)= U21CJ- CJ-p g23C7+ , (23)

where a* = ^ ((7,- ± iOy) are ladder operators. This Hamiltonian anti-commutes with
K = G- and thus mixing preserves the topological charge in ( fTSI ). At = py =
it is independent of the spin rotations up to a global phase of the coupling constants.
Under spin rotation by angle Q all elements in the upper triangular matrix are mul-
tiplied by e'^, while all elements in the lower triangular matrix are multiplied by
e^'^ . This symmetry of triangular matrices does not allow the multiple Fermi point
to split at /7 = 0, as a result the gapless branch of spectrum in Fig. |8]has the cubic
form at low energy, £ -> 0, which corresponds to the topological charge A' = +3:

Topology of quantum vacuum

21

In the low-energy limit the spectrum in the vicinity of the multiple Fermi point
([24t is symmetric under rotations. But in general the spectrum is not symmetric
as demonstrated in Figs. |8] There is only the symmetry with respect to reflection,
{px,Py) —J- {Px,—Py)- The rotational symmetry of spectrum (l24li is an emergent phe-
nomenon, which takes place only in the limit p 0. For trilayer graphene this spec-
trum has been discussed in Ref. ISTI .

4.3 Dirac fermions with quartic and higher order spectrum

3

2

-2

-3

Fig. 9 Spectrum of the Hamiltonian ([25) showing quartic dispersion for the lowest two eigen-
values around the point p = 0. Different colors correspond to different eigenvalues, spectmm is
shown as function of p, at p,, = 0. The spectra have been calculated with equal coupling strengths

gn = gl3 = g22 = gl4 = g24 = gu = g-

In case of four fermionic species, the mixing which does not produce splitting
of the Fermi point is obtained by the same principle as in (l23l i: all matrix elements
aboye the main diagonal contain only (7+ (or (7^):

^(p)

/ c^ P

821 <y^

8i2(y

op

832<y'
842<y'

8n<y 8i4(y^

823<y^ 824<y'^
842(y' (y 1?

\

(25)

Then again under spin rotation by angle 9 all elements in the upper triangular matrix
are multiplied by e'^, while all elements in the lower triangular matrix are multiplied
by e^'^ . As a result the multiple Fermi point with = 4 is preserved giving rise to
the quartic spectrum in vicinity of the Fermi point:

22

G.E. Volovik

For the tetralayer graphene this spectrum was suggested in Ref. ||82l . Again the rota-
tional symmetry emerges new the Fermi point, but spectrum is not symmetric under
permutations. In fact, there is no such 4x4 matrix that would consist of couplings
described by ladder operators and which would be symmetric under permutations.
In general the Fermi point with arbitrary may give rise to the spectrum

£- = l^p'". (27)

Such spectrum emerges in multilayered graphene 1291 l83l . The discussed symmetry
of matrix elements g„,„ extended to 2N x 2N matrix gives ( |27] l with the prefactor

\g\2\\g2-i\---\gN-\M

Violation of this symmetry may lead to splitting of the multiple Fermi point into
elementary Fermi points - Dirac points with A' = 1 and 'relativistic' spectrum

£■2 oc p2.

The effective Hamiltonian describing fermions in the vicinity of multiple Fermi
point is H = + O^pl, see 1291 |3T1 . An example of the effective Hamiltonian

describing the multiple Fermi point with topological charge in 3+1 systems is

ED

H = a,p, + a-pl + a+p^_ . (29)

This Hamiltonian has the spectrum E'^ ~ pi + p^, which has linear dispersion in
one direction and non-linear dispersion in the others.

This is also apphcable to the vacuum of particle physics. The Lorentz symmetry
prohibits both the splitting of the Dirac points and the non-linear non-relativistic
spectrum. Situation changes if the Lorentz symmetry is viewed as an emergent phe-
nomenon, which arises near the Dirac point (Fermi points with A^ = ± 1 ). In this case
both both splitting of Dirac points and formation of non-linear non-relativistic spec-
trum in the vicinity of the multiple Fermi point are possible, an the choice depends
on symmetry. In both cases the mixing of fermions violates the effective Lorentz
symmetry in the low-energy corner. This phenomenon, called the reentrant viola-
tion of special relativity i67J , has been discussed for Nf ^ 3 fermion families in
relation to neutrino oscillations ll58l . Influence of possible discrete flavor symme-
tries on neutrino mixing has been reviewed in Ref. l85l .

5 Flat bands in topological media

Let us extend the results of the Sec. |4.3| to the limit case of infinite number of lay-
ers, N ^ °° 1841 . If the layers are equivalent and interact only via nearest neighbor
couplings, i.e. the nonzero matrix elements are

gn — g23 — ■ --gN-lM = t ,

(30)

Topology of quantum vacuum

23

Py

Px

c

N = 1

Fig. 10 Nodal spiral generates topologically protecetd flat band on the surface. Projection of spiral
on the surface determines boundary of flat band. At each (p., ,Pv) except the boundary of circle
Px + Py = the system represents the ID gapped state (insulator). At each (px,Py) inside the
circle, the insulator is topological being described by non-zero topological invariant and thus
one has gapless edge state. The manifold of these zero-energe edge state inside the circle fonns the
flat band found in Ea.l l32t .

the low-energy spectrum becomes

In the N ^ °° Umit one obtains that in the lowest energy branch all the fermions
within circumference |p| = f have exactly zero energy.

Here we consider the topological origin of the flat band. When the number of layers
N ^ °°, the quasi two dimensional system transforms to the D = 3 system. We find
that the flat band emerges in the first and the last layers, i.e. on the boundaries of
the 3D system. Its formation is accompanied by the simultaneous formation of the
line of node (zeroes of co-dimension 2) in the bulk material. These two objects,
surface flat band and line of zeroes in bulk, are connected via a special kind of the
topological bulk-surface correspondence: projection of the nodal line to the surface
determines the boundary of the flat band (Fig.fTOli.

(31)

£(A^^oc,|p| <f) =0.

(32)

24

G.E. Volovik

5.1 Topological origin of surface fiat band

To understand the topological origin of this branch and its structure let us consider
the spectrum in the continuous limit. The effective Hamiltonian in the 3-dimensional
bulk system which emerges in the limit of infinite number of layers is the following
2x2 matrix

Here t is the magnitude of the hopping matrix element between the layers in Eq.dSOl).
The energy spectrum of the bulk system

E'^ = [px-fcos(fl/5-)]^ + bv + fsin(fl/?,)]2, (34)

has zeroes on line (see Fig.fTOli:

Px = f cos(ap-) , py = —tsm{ap:) , (35)

which forms a spiral, the projection of this spiral on the plane p- = const being the
circle = p^ + Py = t^. This nodal line is topologically protected by the same
topological invariant as in Eq.(fT6]l

N ^ — tr <£ dl a-H-^ViH, (36)
47ti Jc

where the integral is now along the loop C around the nodal line in momentum
space, see Fig. [TO] The winding number around the element of the nodal line is

N=l.

Let us consider now the momentum px as a parameter of the ID system, then
for IpxI 7^ t the system represents the fully gapped system - ID insulators. This
insulator can be described by the same invariant as in Eq.([36]l with the contour of
integration chosen parallel to p-, i.e. along the ID Brillouin zone at fixed px (due to
periodic boundary conditions, the points p^ ^±n/a are equivalent and the contour
of integrations forms the closed loop):

A^(P±) = r-t'" / dp,a,H-Wp_H. (37)

4-711 J-K/a

For |px| < t the ID insulator is topological, since A^(|pxl <t) = 1, while for |px| > t
one has A^i (px ) = and the ID insulator is the trivial band insulator The line |px | =
/ thus marks the topological quantum phase transition between the topological and
non-topological ID insulators.

Topological invariant A^(px) in ( |37] ) determines the property of the surface bound
states of the ID system at each px. Due to the bulk-edge correspondence, the topo-
logical ID insulator must have the surface state with exactly zero energy. Since such
states exist for any parameter within the circle |px | ~ t, one obtains the flat band of
surface states with exactly zero energy, £ ( | Px | <t)=0, which is protected by topol-

Topology of quantum vacuum

25

-1

-0.5

0.5

Px/t

Fig. 11 Formation of the surface flat band. When the number A' of layers increases, the disper-
sionless band evolves from the gapless branch of the spectrum, which has the form E = ±|p^|'^
in the vicinity of multiple Dirac point. The spectrum is shown as a function of for Py = 0. The
curves for A' = 100 and A' = 200 are almost on top of each other. Asymptotically the spectrum
E = ±|pj^|'^ transforms to the dispersionless band within the projection of the nodal line to the
surface.

ogy. This is the origin of the unusual branch of spectrum in Eq.(l32b: it represents the
band of topologically protected surface states with exactly zero energy. Such states
do not exist for parameters |px | > t, for which the ID insulator is non-topological.

The zero energy bound states on the surface of the system can be obtained di-
rectly from the Hamiltonian:

H = Ox{px-tco^{ap,)) + ay{py + t^\n{ap;)) , p^ = —id, , z<0. (38)

We assumed that the system occupies the half-space z < with the boundary at
Z = 0. This Hamiltonian has the bound state with exactly zero energy, ^(px) = 0,
for any |px| < t, with the eigenfunction concentrated near the surface:

The normalizable wave functions with zero energy exist only for px within the circle
IPxI < U i-C- the surface flat band is bounded by the projection of the nodal spiral
onto the surface. Such correspondence between the flat band on the surface and lines
of zeroes in the bulk has been also found in Ref. 1321 for superconductors without
inversion symmetry.

(39)

26

G.E. Volovik

1.5r

N=3

-0.5

-I.5L

-1.5 -1 -0.5 0.5 1 1.5

Px/t

Fig. 12 Fomiation of the nodal line from the evolution of the gapped branch of the spectrum of the
multilayered system, when the number N of layers increases. The spectrum is shown as function of
Px for py = 0. The curves for N = 100 and N = 200 lie almost on top of each other, indicating the
bulk limit. Asymptotically the nodal line p, = fcos(ap.), Py = — fsin(ap;) is formed (two points
on this line are shown, which correspond to Py = 0).

5.2 Dimensional crossover in topological matter: Formation of the
flat band in multilayered system

The discrete model with finite number of layers has been considered in Ref . Il84ll .
It is described by the 2N x 2N Hamiltonian with the nearest neighbor interaction
between the layers in the form:

////px) = (T-p±5,7-f(7+5,j+i -fc7"5,j_i , l<i<N , Px = (/7^,/7v). (40)

In the continuous limit of infinite number of layers ( l40b transforms to (l33T l with
the nodal line in the spectrum, while for finite the spectrum contains the Dirac
point with multiple topological charge equal to A^. Figures (fTTT i - (fl4] | demonstrate
how this crossover from 2D to 3D occurs. When the number A^ of layers increases,
the dispersionless surface band evolves from the gapless branch of the spectrum
E = ijpxl'^. Simultaneously, the gapped branches of the spectrum of the finite-A^
system give rise to the nodal line in bulk. This scenario of formation of the surface
flat band takes place if the symmetry does not allow the splitting of the multiple
Dirac point, and it continuously evolves to the dispersionless spectrum.

If this symmetry is absent, but the symmetry supporting the topological charge
persists, the scenario of the flat band formation is different but still is governed by
topology. Fig. ( fTsT i demonstrates the formation of the flat band in this situation. The
reason, why the flat band emerges is the formation of the topologically protected
nodal line in bulk. The projection of the nodal line on a surface gives the boundary

Topology of quantum vacuum

27

Fig. 13 The lowest energy states for different p., and p,. and arbitrary p- in the bulk limit. The flat
band of surface states is formed in the region p\+p^, <t^. Outside this region only the bulk states
exist. In a simple model considered here flat band comes from the degenerate Dirac point with
nonlinear dispersion. However, this is not necessary condition: the flat band emerges whenever the
nodal line appears in the bulk.

Py/i

Fig. 14 Result of the transformation of the gapped states in the process of dimensional crossover.
They form the nodal line in bulk p, =?cos(op.), p,, = — f sin(ap-) whose projection to the (p, , p, )-
plane is shown.

of the flat band emerging on this surface. This is the realization of the bulk-surface
correspondence in systems with the nodal lines in bulk. The other examples can be
found in Ref. US).

28

G.E. Volovik

N-ino

N = :200

-1.5

-1

-0.5

0.5

1

1.5

P. It

Fig. 15 Evolution of the spectrum in a different model, in which the multiple Dirac point is absent,
but the suiface flat band is formed together with the formation of the nodal line in bulk

6 Anisotropic scaling and Hof ava gravity

We know that in the vicinity of the Weyl point with elementary topological charge
A^ = +l orA^ = — 1 the quantum electrodynamics and gravity emerge as effective
fields. Here we discuss what kind quantum electrodynamics emerges near the node
with higher A^. This is the quantum electrodynamics with anisotropic scaling, in
which as distinct from isotropic relativistic QED the electric and magnetic fields
obey different scaling laws. The analogous anisotropic scaling has been discussed
for quantum gravity at short distances.

6.1 Effective theory near the degenerate Dirac point

Let us consider again the 2+1 system with the multiple Dirac point and write the
following effective Hamiltonian:

J^N = ((ei +/e2) • (p-eA))l"l + ((ei - /e2) • {v-eP^)P .

For a single layer (A^ = 1) the Hamiltonian ( |4TI ) is reduced to the conventional Weyl-
Dirac Hamiltonian for massless particles in 2+1 dimension;

=C7-'ei-(p-eA) + (7-^'e2-(p-eA)=4c7"(/7;-eA,) , fl = (l,2). (42)

Here the vectors ei and 62 play the role of zweibein (in the ground state they are
mutually orthogonal). The vector A is either the vector potential of the effective
electromagnetic field which comes from the shifts of the node, or the real electro-

(41)

Topology of quantum vacuum

29

magnetic field, as it takes place for electrons in graphene. The whole dreibein e„
with a ~ (1,2,3) and = (0, 1,2) emerges for the Green's function and gives the
effective metric 2+1 metric as a secondary object:

g^' = il"'e^el , (43)

For general |A^| situation is somewhat different. While the action for the relativis-
tic fermions is invariant under rescaling r = br', t = bt' , the action for the fermions
living in the vicinity of the multiple Dirac point is invariant under anisotropic rescal-
ing r = br' , t = Zjl^lf'. The anisotropic scaling is in the basis of the Hofava gravity,
which is described by the space components of metric are separated from the time
component and have different scaling laws [62 63 MllSl- The square of Hamilto-
nian (4T[ gives the space metric in terms of zweibein:

jr^2^4=(g'^(p,-M,)(p,-M,))'"' , g'^ =e\ei+e^ei. (44)

6.2 Effective electromagnetic action

If the effective action for fields Ca and is obtained by integration over fermions in
the vicinity of the multiple Dirac point, this bosonic action (actually the terms in ac-
tion which mostly come from these fermions) inherits the corresponding conformal
symmetry of the massless fermions. For anisotropic scaling the conformal invari-
ance means invariance under g''^ — > b^g''' and g^^ — > b^^^^g^^; y/—g
(where D is space dimension); while g^' is not considered.

6.2.1 Single layered graphene and relativistic fields

For D y^2 the spectrum of multiple Fermi point becomes more complicated, and in
general is not isotropic, see Eq.(l29]l. For general D the spectrum is isotropic only
for |A^| = 1, where one obtains effective relativistic massless D+l quantum electro-
dynamics, which is Lorentz invariant. This implies the following nonlinear action

Se,„{\N\ 1,£>) = j d'^xdt [B^ _£2j ^ ^ (-45-)

For D = 3 the action is proportional to (Z?^ — E^)\n{B^ — E^), and is imaginary at
< giving rise to Schwinger pair production in massless quantum electrody-
namics. The similar imaginary action tales place for < E^ for D ^3. For exam-
ple, for a single layer graphene (D = 2, |A^| = 1) reproduces the relativistic 2+1 QED
which gives rise to Lagrangian (B^ — fi^)^/^ fST! with the running coupling constant
1/a = {3/ 2)/ 871^. The action is imaginary at B^ < E^ which corresponds to
Schwinger pair production with the rate E^^^ at B = 0.

30

G.E. Volovik

6.2.2 Bilayer graphene

For bilayered graphene, assuming the quadratic dispersion |A^| = D = 2, the expected
conformal invariant Heisenberg-Euler action for the constant in space and time elec-
tromagnetic field, which is obtained by the integration over the 2+ 1 fermions with
quadratic dispersion, is the function of the scale invariant combination jj. Il88l[89l :

£■2

d^xdt B^gipL) , ^ = (46)

The asymptotical behavior in two limit cases, ^(/i 0) ~ const and — > oo) ~
^2/3^ gives the effective actions for the constant in space and time magnetic and
electric fields:

Sb

I J dhdt , SE^{b + ic) j d^xdtE^I^ . (47)

The parameter a is the logarithmic coupling constant; the parameter b describes the
vacuum electric polarization; and the parameter c describes the instability of the
vacuum with respect to the Schwinger pair production in the electric field, which
leads to the imaginary part of the action. The action also contains the linear non-
local term

Snon-locaim ^ D ^ 2) = J d" xdt ^g^' g'" Fak—^FQ„ , (48)

which corresponds to the polarization operator

7100--== , \N\=D = 2. (49)
V^* - or-

6.2.3 D=2 systems with nodes with topological charge A'

In general case of a 2D system with A^-th order touching point in spectrum (a kind of
multilayered graphene) the Heisenberg-Euler action contains among the other terms
the following nonUnear terms in the actions for magnetic and electric fields

Sb{N,D^2)^ Jd^xdtB^ , Se{N,D = 2)^ j d^xdt{-E^Y^(^+W) , (50)

where imaginary part of the action is responsible for the pair production in electric
field |i88„90J, while the Hnear action of the type ( |48] | corresponds to the polarization
operator

Topology of quantum vacuum

6.2.4 Effective action for gravity

The expected action for gravitational field is

31

Sgrilv — d xdt^/—g

(52)

where Ki,K2,... are dimensionless quantities, and the other terms of that type are
implied.

7 Fully gapped topological media

Examples of the fully gapped topological media are the so-called topological band
insulators in crystals [3]. Examples are Bi2Se3, Bi2Te3 and Sb2Te3 compounds
which are predicted to be 3+1 topological insulators |91j. But the first discussion
of the 3+1 topological insulators can be found in Refs. 1921 |93l . The main feature
of such materials is that they are insulators in bulk, where electron spectrum has
a gap, but there are 2+1 gapless edge states of electrons on the surface or at the
interface between topologically different bulk states as discussed in Ref. ll93l . The
similar properties are shared by the fully gapped 3D topological superfluids and su-
perconductors. The spin triplet p-wave superfluid ^He-B represents the fully gapped
superfluid with nontrivial topology. It has 2+1 gapless quasiparticles living at inter-
faces between vacua with different values of the topological invariant describing
the bulk states of ^He-B lll4ll94ll . The quantum vacuum of Standard Model below
the electroweak transition, i.e. in its massive phase, also shares the properties of the
topological insulators and gapped topological superfluids and is actually the rela-
tivistic counterpart of ■'He-B 1231 .

Examples of the 2+1 topological fully gapped systems are provided by the films
of superfluid ^He-A with broken time reversal symmetry 12111951 and by the planar
phase which is time reversal invariant l2ni95]| . The topological invariants for 2+1
vacua give rise to quantization of the Hall and spin-Hall conducticity in these films
in the absence of external magnetic field (the so-called intrinsic qauntum and spin-
quantum Hall effects) 12111961 .

7.1 2+1 fully gapped vacua

7.1.1 ^He-A film: 2+1 chiral superfluid

The gapped (nodeless) ground states (vacua) in 2+1 systems are characterized by
the invariant obtained by dimensional reduction from the topological invariant de-
scribing the nodes of co-dimension 3. The invariant for the Fermi point in (|2|is the
integral over the 3D surface a around the singularity. The invariant describing the

32

G.E. Volovik

Fig. 16 Skyrmion in p-space with momentum space topological charge A' = — 1 in (55). It de-
scribes topologically non-trivial vacua in 24-1 systems with a fully gapped non-singular Green's
function. Vacua with nonzero N have topologically protected gapless edge states. The nonzero
topological charge leads also to quantization of Hall and spin Hall conductance.

2D insulator or nodeless supefluid is the integral over the whole (2+l)-dimensional
momentum-frequency space {px-Py, co):

24n-

d^pdco Gdp. G-^GdpG-^ Gdp^ G '

(53)

This equation ( l53T l was introduced in relativistic 2-1-1 theories ifTTl [121 [131 and
for the film of ^He-A in condensed matter |[T6l [2TI . where it was inspired by the
dimensional reduction from the Fermi point, see 1951 . In simple case of the 2x2
matrix, the Green's function can be expressed in terms of the three-dimensional
vector d(/7^,/7y),

G-\(0,px.p,)^i(0-H , H^x-d{p,„py), (54)

where H is the Bogoliubov - de Gennes Hamiltonian for fermions in p-wave super-
fluids, and T are the Pauli matrices. Example of the d-vector configuration, which
corresponds to the topologically nontrivial vacuum is presented in Fig. [16] This is
the momemtum-space analog of the topological object in real space - skyrmion. In
real space, skyrmions are described by the relative homotopy groups 1971 ; they have
been investigated in detail both theoretically and experimentally in the A phase of
^He, see Sec. 16.2 in 11] and the review paper 1981 .

For the Green's function in i54[ the winding number of the momentum-space
skyrmion in Eq.([53]l is reduced to lT6l

N^T- d^pd-i^x — ] , (55)

47t J \dpx dpy

where d = d/|d| is unit vector. For a single layer of the ^He-A fibn and for one spin
projection, the simplified Bogoliubov - de Gennes Hamiltonian has the form:

Topology of quantum vacuum

33

N = 4

N = 2

gap

N = Ql\

gap

gap

a

— ►

trivial topological superfluids ^1™ thickness

superfluid

Fig. 17 Dependence of the topological invariant i53\ on the thickness of ^^He-A film. The even
values of N result from the spin degeneracy. At the topological phase transitions between the states
with different A', the gap in the spectrum of fermions is nullified. The interface between the vacua
with different topological charges contains gapless fermions - edge states. The number of the
gapless fermions is related to the difference of the topological charges by index theorem.

H = T- d(p) = iCO + T3 y-^^^ - M j + + TlPy , (56)

For ;Li > the topological charge in Eq.dSSTi or in Eq.(l53]l is = 1, and the d(p)
field forms the skyrmion in momentum space of the type shown in Fig. [16] For
jj. < the topological charge is trivial, N = 0. That is why at /i = there is a
topological quantum phase transition between the topological superfluid at /i >
and non-topological superfluid at /i < 1951 .

In general case of multilayered ^He-A, topological charge A^ may take any integer
value of group Z. This charge determines quantization of Hall and spin-Hall conduc-
tance and the quantum statistics of the topological objects - real-space skyrmions
|[l6l|2ll|96l|95l. For A^ = 4;t + 1 and A^ = 4yt + 3, skyrmion is anyon; for A^ = 4yt + 2
it is fermion; and for A^ = 4A: it is boson ||95|| . This demonstrates the importance
of the Z2 and Z4 subgroups of the group Z in classification of topological matter;
and also provides an example of the interplay of momentum-space and real-space
topologies.

7.1.2 Planar phase: time reversal invariant gapped vacuum

In case when some symmetry is present, additional invariants appear, which corre-
spond to dimensional reduction of topological invariant supported by symmetry Nk
in®:

-1

j cf-pd(0 KGdpfi- ' Gdpj G ' Gdp^G-

(57)

34

G.E. Volovik

where the matrix K commutes with the Green's function matrix. Example of the
symmetric 2+1 gapped state with Nk is the film of the planar phase of superfluid
^He 1211 |95l . In the single layer case, the simplest expression for the Green's func-
tion is

where H is the Bogoliubov - de Gennes Hamiltonian for fermions in this spin-triplet
p-wave superfluid. The symmetry operator K, which supports the topological in-
variant \sK ~ 13(7-: it commutes with the Green's function. The planar state is time
reversal invariant. It has trivial conventional topological charge N ^0 and non-zero
symmetry protected charge Nk ~ 2. For the general case of the quasi 2D film with
multiple layers of the planar phase, the invariant Nk belongs to the group Z. The
magnetic solid state analog of the planar phase is the 2D time reversal invariant
topological insulator, which experiences the quantum spin Hall effect without exter-
nal magnetic field ||3l .

8 Relativistic quantum vacuum and superfluid ^He-B

Let us now turn to the class of 3+1 fully gapped systems, which is represented by
Standard Model in its massive phase and superfluid ^He-B.

In the broken symmetry phase of Standard Model below the electroweak transi-
tion, symmetry does not support the topological invariant responsible for the nodes
in spectrum. In this phase there is no mass protection by topology and thus all the
fermions become massive, i.e. Standard Model vacuum becomes the fully gapped
insulator

In quantum liquids, the fully gapped three-dimensional system, which is simi-
lar to the vacuum of massive Standard Model, is represented by another phase of
superfluid ■'He - the ^He-B. This phase has time reversal symmetry and nontrivial
topology supported by symmetry, which gives rise to the 2D gapless quasiparticles
living at interfaces between vacua with different values of the topological invariant
or on the surface of ^He-B |[T4l |94l |99l [TOOl .

8.1 Superfluid ^He-B

^He-B belongs to the same topological class as the vacuum of Standard Model in
its present insulating phase ||23| . The topological classes of the ^He-B states can be
represented by the following simplified Green's function and the Bogoliubov - de
Gennes Hamiltonian:

(58)

Topology of quantum vacuum

35

G (co,p) = i(o-H , // = T3 M I + Tic^a p.

' 2m

(59)

In the limit 1/m = this model ^He-B transforms to the vacuum of massive rela-
tivistic Dirac particles with speed of light c = and mass parameter M = —ji.

In the fully gapped systems, the Green's function has no singularities in the whole
4-dimensional space (ft),p). That is why we are able to use the Green's function
at (B = 0, which corresponds to the effective Hamiltonian, i/eff(p) = ^G^'(0,p).
The topological invariant relevant for ^He-B and for quantum vacuum with massive
Dirac fermions is:

247r2

tr

(B=0

d^'p K Gdpfi- ' Gdp G- ' Gdpfi

— 1

(60)

with matrix = T2 which anti-commutes with the Green's function at = 0. In
^He-B, the T2 symmetry is combination of time reversal and particle-hole symme-
tries; for Standard Model the matrix T2 = 757". Note that at CO = the symmetry
of the Green's function is enhanced, and thus there are more matrices K, which
commute or anti-commute with the Green's function, than at O) 7^ 0.

1/m

non-topological superfluid

topological superfluid 'hc-B

N" =0

N" =+2

i

n'' =-\

N*^ =+\

y

Dime

K

Dirac

N" =-1

N" =0

P topological superfluid

non-topological superfluid

Fig. 18 Phase diagram of topological states of ■'He-B in Eq. <59| l in the plane [p., 1/m). States on
the line 1 /m = correspond to the Dirac vacua, which Hamiltonian is non-compact. Topological
charge of the Dirac fermions is intemiediate between charges of compact ^He-B states. The line
l/m — separates the states with different asymptotic behavior of the Green's function at infinity:
G^' ((B = 0,p) — > ±tifp- /2m. The line = marks topological quantum phase transition, which
occurs between the weak coupling ■'He-B (with > 0, ;;j > and topological charge Nk = 2)
and the strong coupling ■'He-B (with /i < 0, m > and Nk = 0). This transition is topologically
equivalent to quantum phase transition between Dirac vacua with opposite mass pai'ameter M =
it|/i|, which occurs when /i crosses zero along the line 1/m = 0. The interface which separates
two states contains single Majorana fermion in case of ^He-B, and single chiral fermion in case
of relativistic quantum fields. Difference in the nature of the fermions is that in Fermi superfluids
and in superconductors the components of the Bogoliubov-Nambu spinor ai'e related by complex
conjugation. This reduces the number of degrees of freedom compared to Dirac case.

36

G.E. Volovik

Fig. [18] shows the phase diagram of topological states of ^He-B in the plane
{jj., 1 /m). The line 1/m = corresponds to the Dirac vacuum of massive fermions,
whose topological charge is determined by the sign of mass parameter M = —jj.:

Nk = sign{M). (61)

The real superfluid ^He-B lives in the corner of the phase diagram jj.>0, m > 0,
fi ^ mc\, which also corresponds to the limit of the weakly interacting gas of ^He
atoms, where the superfluid state is described by Bardeen-Cooper- Schrieffer (BCS)
theory. However, in the ultracold Fermi gases with triplet pairing the strong coupling
limit is possible near the Feshbach resonance 11011 . When pL crosses zero the topo-
logical quantum phase transition occurs, at which the topological charge AOs: changes
from Nk = 2 to Nk = 0. The latter regime with trivial topology also includes the
Bose-Eistein condensate (BEC) of two-atomic molecules. In other words, the BCS-
BEC crossover in this system is always accompanied by the topological quantum
phase transition, at which the topological invariant changes.

There is an important difference between ^^He-B and Dirac vacuum. The space of
the Green's function of free Dirac fermions is non-compact: G has different asymp-
totes at IpI oo for different directions of momentum p. As a result, the topologi-
cal charge of the interacting Dirac fermions depends on the regularization at large
momentum. ^He-B can serve as regularization of the Dirac vacuum, which can be
made in the Lorentz invariant way 1231 . One can see from Fig. [18] that the topo-
logical charge of free Dirac vacuum has intermediate value between the charges of
the ^He-B vacua with compact Green's function. On the marginal behavior of free
Dirac fermions see Refs. |[T5] [19] [111991 .

The vertical axis separates the states with the same asymptote of the Green's
function at infinity. The abrupt change of the topological charge across the line,
ANk = 2, with fixed asymptote shows that one cannot cross the transition line adi-
abatically. This means that all the intermediate states on the line of this QPT are
necessarily gapless. For the intermediate state between the free Dirac vacua with
opposite mass parameter M this is well known. But this is applicable to the gen-
eral case with or without relativistic invariance: the gaplessness is protected by the
difference of topological invariants on two sides of transition.

8.2 From superfluid relativistic medium to ^He-B

We have already seen, that the Dirac vacuum of massive Standard Model particles
has nontrivial topology. As a result the domain wall separating vacua with opposite
signs of the mass parameter M contains fermion zero modes II102I . Other examples
of the topologically nontrivial states in relativistic theories can be provided by dense
quark matter, where chiral and color superconductivity is possible. The topological
properties of such fermionic systems have been recently discussed in Ref. 110311 . In
particular, in some range of parameters the isotropic triplet relativistic superconduc-

Topology of quantum vacuum

37

tor is topological and has the fermion zero modes both at the boundary and in the
vortex core. On the other hand, there is a range of parameters, where this triplet su-
perconductor is reduced to the non-relativistic superfluid ^He-B 11041 . That is why
the analysis in Ref. 11031 is applicable to ^He-B and becomes particularly useful
when the fermions living in the vortex core are discussed.

In relativistic superconductor or superfluid with the isotropic pairing - such as
color superconductor in quark matter - the fermionic spectrum is determined by
Hamiltonian

H = T3{ca-p + pM- iir) + TiA , (62)
for spin singlet pairing, and by Hamiltonian

// = T3(ca-p + j3M-jLiff) + 75Ti4, (63)

for spin triplet pairing II1041 11031 . Here a', jS and 75 are Dirac matrices, which in
standard representation are

M is the rest energy of fermions; is their relativistic chemical potential as dis-
tinct from the non-relativistic chemical potential jj.; T„ are matrices in Bogoliubov-
Nambu space; and A is the gap parameter.

In non-relativistic limit the low-energy Hamiltonian is obtained by standard pro-
cedure, see e.g. 11051 . The non-relativistic limit is determined by the conditions

cp < M (65)

and

\M-^J^l + A^\<^M. (66)

Under these conditions the Hamiltonian ( l62b reduces to the Bogoliubov - de Gennes
(BdG) Hamiltonian for fermions in spin-singlet i-wave superconductors, while ( l63l l
transforms to the BdG Hamiltonian relevant for fermions in isotropic spin-triplet
/7-wave superfluid ^He-B in Eq.

, ' \ R M „ A

where the nonrelativistic chemical potential jj. = + A^ — M.
The Dirac-BdG system in Eq.(l63]l has the following spectrum

e = ±]jM^+c^p^ + A^ + lil±2^M^{nl + A-)+n^c^-p^. (68)

This spectrum is plotted in Fig(T9] Depending on the value of the parameters Hr, A,
M the spectral branches have different configurations.

38

G.E. Volovik

Fig. 19 Plot of the spectrum of relativistic Hamiltonian <63b for two generic cases: (a) \i1r\ > jl^
when the minimum in the energy spectrum is away from the origin and (b) < fl^ when the
minimum in the energy spectrum is at p = 0. At > /i^ there is a soft quantum phase transition
between these two vacua. This transition is not topological, and thus the gap in the energy spectrum
does not close at the transition. The gap closes at the topological transition occurring at fi^ =
— as shown in plot (c).

There is a soft quantum phase transition, at which the position of the minimum of
energy E(p) shifts from the origin p = 0, and the energy profile forms the Mexican
hat in momentum space. This momentum-space analog of the Higgs transition ||2l
occurs when the relativistic chemical potential exceeds the critical value

^*^=[^ + \l^+^^') ' ■ ^^^^

Figures [19] (a) and[T9](b) demonstrate two generic cases: l/i^j > /i^ when there are
extremums of function e(p) at /9 7^ and \plr\ < jJ-g when all extremums are at
the point p = 0. The formation of the Mexican hat at = /i^ is an example of
non-topological quantum phase transition, which occurs without change of the topo-
logical invariants and thus is not accompanied by the gapless intermediate state. Let
us turn to the topological quantum phase transitions, at which the p-space topolog-
ical invariant changes and the gap closes at the transition point as is shown in Fig.
[11(c).

8.3 Topology of relativistic medium and ^He-B

Fig. [20] shows the phase diagram of the vacuum states of relativistic triplet super-
conductors. Different vacuum states are characterized by different values of the
topological invariant Nk in Eq.(i60]l, where the Green's function matrix at zero fre-

Topology of quantum vacuum

39

quency G^' (o) = 0,p) is equivalent to effective Hamiltonian. For ^He-B in Eq.(l67]i
and for triplet relativistic superconductor in Eq. (|63] | the relevant matrix = T2,
which anti-commutes with the Hamiltonian. The vacuum states with different Nk
cannot be adiabatically connected, and thus at the phase transition lines the states
are gapless. The circle /xj + 4^ = is an example of the line of topological quan-
tum phase transition. In non-relativistic limit this corresponds to the line /i = in
Figim The states inside the circle + = are topologically trivial, while the
states outside this circle represent topological superconductivity II103II . The vacuum
states with /ij + 4^ > and /ij +4^ < cannot be adiabatically connected
which leads to the gap closing in Fig[T9](c). Discontinuity in the topological charge
across the transition induces discontinuity in the energy of the ground state across
the transition. For example, for the 2+1 px + ipy superfluid/superconductor the quan-
tum phase transition is of third order, meaning that the third-order derivative of the
ground state energy is discontinuous II106I .

Fig. 20 Phase diagram of ground states of relativistic triplet superfluid in Eg. <63t in the plane
{Hr,A). Topological quantum phase transitions are marked by thick lines. The states inside the
circle pL^ + A- = M- are topologically trivial. The states outside this circle represent topological
superconductors. The states on the lines of topological quantum phase transition are gapless.

9 Fermions in the core of strings in topological materials

In relativistic theories there is an index theorem which relates the number of fermion
zero modes localized on a vortex with the vortex winding number |10?1. We know
that the Dirac vacuum considered in Ref. 11071 has nonzero topological charge.
This suggests that the existence of zero energy states in the core is sensitive not only
to the real-space topological charge of a vortex, but also to the momentum-space
topological charge of the quantum vacuum in which the vortex exists, and if so the
index theorem should be extended to vortices in any fully gapped systems, including
the non-relativistic superfluid ^He-B. Here we discuss this issue of the connection

40

G.E. Volovik

between the topological charge of the vacuum and existence of the fermion zero
modes on the topological objects in this vacuum. Our examples demonstrate that
during the topological quantum phase transition from the topologically nontrivial
vacuum to the trivial one, the fermion zero modes on vortices (or other objects)
disappear.

9.1 Vortices in ^He-B and relativistic strings

Since in this problem both the momentum-space topology of bulk state and the
real-space topology of the vortex or other topological defects are involved, the com-
bined topology of the Green's function in the coordinate-momentum space (©, p, r)
ll38l 11081 [T] 11091 II 101 seems to be relevant. However, though the bulk- vortex cor-
respondence does evidently exist, the explicit index theorem which relates the exis-
tence of the fermion zero modes to the topological charge of the bulk state and the
vortex winding number is still missing. The existing index theorems are applicable
only to particular cases, see e.g. II1081 11031 II 101 [TTTI . There is also a special in-
dex theorem for superconductors/superfluids with a small gap A pL. Spectrum of
fermions in these superconductors has branches which cross zero energy as a func-
tion discrete quantum number - angular momentum L 11 121 . The index theorem
relates the number of such branches with the vortex winding number 11 131 . Here
we are interested in the true fermion zero modes - the branches of spectrum £(/?-),
which cross zero as function of momentum along the vortex line.

Example, which demonstrates that the connection between the topological charge
Nk and the existence of Majorana fermions - fermion zero modes on vortices - is
in Fig. |2T] For ^He-B, which lives in the BCS range of parameters where Nk ^ 0,
the gapless fermions in the core have been found in Ref. [ 114|. On the other hand,
in the strong coupling limit the ^He-B transforms to the Bose-Einistein condensate
(BEC) of molecules. The latter does not contain fermionic excitations, and thus one
should not expect the existence of gapless fermions in the vortex core. Thus one ex-
pects that somewhere in the region of crossover between the BCS-like and BEC-like
regime the spectrum of fermions localized on vortices must be reconstructed. On the
other hand, the topological reconstruction of the fermionic spectrum in the vortex
core cannot occur during the adiabatic deformation. The discontinuous deformation
of the spectrum in the core is only possible during the topological quantum phase
transition in bulk: at such transition the intermediate bulk gapless state is crossed
which destroys the adiabaticity. This is just what happens: in the BEC limit, the
chemical potential /i is negative and the topological charge of the vacuum state is
trivial, Nk = 0. The reconstruction of the bulk spectrum at the topological quan-
tum phase transition occuring at /i = 0, see Fig. [18] triggers reconstruction of the
spectrum of fermion zero modes in the core. At < the topological charge Nk
nullifies and simultaneously the gap in the spectrum of core fermions arises, see
Fig. |2T| The similar situation, when the quantum phase transition in bulk produces
leads to appearance or disappearance of fermion zero modes, has been discussed in

Fig. 21 Schematic illustration of spectrum of the fermionic bound states in the core of the
most symmetric vortex with « = 1, the so-called o- vortex 1981 . in fully gapped spin triplet su-
perfluid/superconductor of ■'He-B type, (left): Spectrum of bound state in the ■'Hc-B, vortex which
corresponds to the weak coupling BCS regime with non-zero topological charge Nk = 2 11141 . L is
the azimuthal quantum number of fermions in the vortex core. There are two fermion zero modes,
which cross zero energy in the opposite directions, (right): The same vortex but in the topologically
trivial BEC regime, = 0, does not have fermion zero modes. The spectrum of bound states is
fully gapped. Fermion zero modes disappear at the topological quantum phase transition, which
occurs in bulk liquid at ;U = 0. Similar situation may take place for strings in color superconductors
in quark matter 11031 .

Ref. Ill 151 for the other type of p-wave vortices; in Ref. II103I for strings in color
superconductors in quark matter; and in Refs. 01 16111 17l for Majorana fermions on
the edges of quantum wire (the review on Majorana fermions in superconductors
can be found in Ref. 11 181 ).

Another example is provided by the fermions on relativistic vortices in Dirac
vacuum discussed in Ref. II107II . The Dirac vacuum has the nonzero topological in-
variant, Nk = ±1, see Fig. [18] This is consistent with the existence of the fermion
zero modes on vortices, found in Ref. 11071 . The index theorem for fermion zero
modes on these vortices can be derived using the topology in combined coordi-
nate and momentum space. The number of fermion zero modes on a vortex A'zm
can be expressed via the 5-form topological invariant in terms of Green's function
G(£0,p,r) EOlinSl

42

G.E. Volovik

(70)

The space integral is along the closed contour C around the vortex line. For the
vortex in Dirac vacuum, equation dTOb reproduces the index theorem discussed in
Ref. II107II : the algebraic number of fermion zero modes equals the vortex winding
number A^zm = «•

For vortex in ^He-B one obtains Nzm ~ 0. This is, however, consistent with
Fig. |2T] two branches of zero modes have opposite signs of velocity — dE / d p,.
Though the algebraic sum of zero modes, A^zm = 1 — 1=0, due to the special sym-
metry of the vortex configuration the two branches of fermion zero modes do not
cancel each and the gap is not formed. To resolve the fermion zero modes in the sys-
tems with symmetry, the index theorem for the zero modes must be supplemented
by symmetry consideration.

The 5-form topological invariant similar to Eq.dTOli has been discussed also in
Ifn il20l . In particular it is responsible for the topological stability of the 3h-1 chiral
fermions emerging in the core of the domain wall separating topologically differ-
ent vacua in A+\ systems (see Sec. 22.2.4 in yj). The topological invariant for
the general 2n + 1 insulating relativistic vacua and the bound chiral fermion zero
modes emerging there have been considered in 01211 11221 l60l . Application of the
5-form topological invariant to the states in lattice chromodynamics can be found in

9.2 Flat band in a vortex core: analog of Dirac string terminating
on monopole

E(Pz)

continuous spectrum

Pz

Fig. 22 Schematic illustration of the spectrum of bound states E(p~) in the vortex core. The
branches of bound states terminate at points where their spectrum merges with the continuous
spectrum in the bulk. The flat band terminates at points where the spectrum has zeroes in the bulk,
i.e. when it merges with Weyl points. It is the p-space analog of a Dirac string terminating on a
monopole, another analog is given by the Fermi arc in Fig. 1 bottom right .

Topology of quantum vacuum

43

The topological bulk-vortex correspondence exists also for vortices in gapless
vacua. The topological protection of fermion zero modes is provided by the non-
trivial topology of three-dimensional Weyl points in the bulk. This bulk-vortex cor-
respondence II 1001 is illustrated in Fig. |22l In bulk there is a pair of Weyl points
with opposite topological charges = ±1 in Eq.©. The projections of these Weyl
points on the direction of the vortex line determine the boundaries of the region
where the spectrum of fermions bound to the vortex core is exactly zero, E{py ) ~ 0,
for all within this region. Such flat band was first obtained in Ref. 11241 for the
noninteracting model. However, due to the topological protection, it is not destroyed
by interactions. The spectrum of bound states in a singly quantized vortex in ^He-A
is shown in Fig. |22] The ID flat band terminates at points where the spectrum of
bound state merges with zeroes in the bulk, i.e. with Weyl points.

10 Discussion

The last decades demonstrated that topology becomes a very important tool in
physics. Topology in momentum space is the main characteristics of the ground
states of a system at zero temperature (T = 0), in other words it is the characteristics
of quantum vacua. The gaplessness of fermions in bulk, on the surface or inside the
vortex core is protected by topology, and thus is not sensitive to the details of the
microscopic physics (atomic or trans-Planckian). Irrespective of the deformation of
the parameters of the microscopic theory, the value of the gap (mass) in the energy
spectrum of these fermions remains strictly zero. This solves the main hierarchy
problem in particle physics: for fermionic vacua with Fermi points the masses of
elementary particles are naturally small.

The vacua, which have nontrivial topology in momentum space, are called the
topological matter, and the quantum vacuum of Standard Model is the representative
of the topological matter alongside with topological superfluids and superconduc-
tors, topological insulators and semi-metals, etc. There is a number of of topological
invariants in momentum space of different dimensions. They determine universality
classes of the topological matter and the type of the effective theory which emerges
at low energy and low temperature. In many cases they also give rise to emergent
symmetries, including the effective Lorentz invariance and probably all the sym-
metries of Standard Model, and emergent phenomena such as gauge and gravita-
tional fields. The symmetry appears to be the secondary factor, which emerges in
the low-energy corner due to topology, and it is possible that it is topology of the
quantum vacuum, which is responsible for the properties of the fermionic matter in
the present low-energy Universe.

The topological invariants in extended momentum and coordinate space deter-
mine the bulk-surface and bulk-vortex correspondence. They connect the momen-
tum space topology in bulk with the real space. These invariants determine the
fermion zero modes living on the surface of a system or in the core of topologi-
cal defects (vortices, strings, domain walls, solitons, hedgehogs, etc.).

44

G.E. Volovik

In respect to gravity, the momentum space topology gives some lessons. First,
in the effective gravity emerging at low energy, the collective variables represent
the tetrad field and spin connections. In this approach the metric field emerges as
the composite object of tetrad field, and thus the Einstein-Cartan-Sciama-Kibble
theory with torsion field is more relevant for the description of gravity (see also
11125111261 [1271 (1281 ).

Second, the topology suggests several scenarios of Lorentz invariance violation
governed by topology. Among them the splitting of Fermi point and development of
the Dirac points with quadratic and cubic spectrum. The latter leads to the natural
emergence of the Hof ava-Lifshitz gravity.

Acknowledgements This work is supported in part by the Academy of Finland and its COE pro-
gram 2006-2011

References

1. G.E. Volovik, The Universe in a Helium Droplet, Clarendon Press, Oxford (2003).

2. G.E. Volovik, Quantum phase transitions from topology in momentum space, in: "Quan-
tum Analogues: From Phase Transitions to Black Holes and Cosmology", eds. W.G.
Unruh and R. Schiitzhold, Springer Lecture Notes in Physics 718 (2007), pp. 31-73;
cond-mat/0601372

3. M.Z. Hasan and C.L. Kane, Topological Insulators, Rev. Mod. Phys. 82, 3045-3067 (2010).

4. Xiao-Liang Qi and Shou-Cheng Zhang, Topological insulators and superconductors. Rev.
Mod. Phys. 83, 1057-1110(2011).

5. H.B. Nielsen, M. Ninomiya: Absence of neutrinos on a lattice. I - Proof by homotopy theory,
Nucl. Phys. B 185, 20 (1981); Absence of neutrinos on a lattice. II - Intuitive homotopy
proof, Nucl. Phys. B 193, 173 (1981).

6. G.E. Volovik and VP Mineev, Current in superfluid Fermi liquids and the vortex core struc-
ture, lETP 56, 579-586 (1982).

7. H.B. Nielsen and M. Ninomiya, The Adler-Bell-Jackiw anomaly and Weyl fermions in a
crystal, Phys. Lett. 130 B, 389-396 (1983).

8. J.E. Avron, R. Seller and B. Simon, Homotopy and quantization in condensed matter
physics, Phys. Rev. Lett. 51, 51-53 (1983).

9. G.W. Semenoff, Condensed-matter simulation of a three-dimensional anomaly, Phys. Rev.
Lett. 53, 2449-2452 (1984).

10. Qian Niu, D. I. Thouless, and Yong-Shi Wu, Quantized Hall conductance as a topological
invariant, Phys. Rev. B 31, 3372-3377 (1985).

11. H. So, Induced topological invariants by lattice fermions in odd dimensions. Prog. Theor
Phys. 74, 585-593 (1985).

12. K. Ishikawa and T. Matsuyama, Magnetic field induced multi component QED in three-
dimensions and quantum Hall effect, Z. Phys. C 33, 41^5 (1986).

13. K. Ishikawa and T. Matsuyama, A microscopic theory of the quantum Hall effect, Nucl.
Phys. B 280, 523-548 (1987).

14. M.M. Salomaa and G.E. Volovik, Cosmiclike domain walls in superfluid ■'He-B: Instantons
and diabolical points in (k,r) space, Phys. Rev. B 37, 9298-9311 (1988).

15. F.D.M. Haldane, Model for a quantum Hall effect without Landau levels: Condensed-matter
realization of the "Parity Anomaly", Phys. Rev. Lett. 61, 2015-2018 (1988).

16. G.E. Volovik, Analog of quantum Hall effect in superfluid ^He film, JETP 67, 1804-181 1
(1988).

Topology of quantum vacuum

45

17. V.M. Yakovenko, Spin, statistics and charge of solitons in (2+l)-dimensional theories, Fizika
(Zagreb) 21, suppl. 3, 231 (1989); arXiv:cond-mat/9703195

18. P. Hofava, Stability of Fermi surfaces and AT-theory, Phys. Rev. Lett. 95, 016405 (2005).

19. A.P. Schnyder, S. Ryu, A. Furusaki and A.W.W. Ludwig, Classification of topological in-
sulators and superconductors in three spatial dimensions, Phys. Rev. B 78, 195125 (2008);
A.P. Schnyder, S. Ryu, A. Fumsaki and A.W.W. Ludwig, Classification of topological insu-
lators and superconductors, AIP Conf. Proc. 1134, 10 (2009); arXiv; 0905. 2029

20. A. Kitaev, Periodic table for topological insulators and superconductors, AIP Conference
Proceedings, Volume 1134, pp. 22-30 (2009); arXiv:0901.2686

21. G.E. Volovik and V.M. Yakovenko, Fractional charge, spin and statistics of solitons in su-
perfluid ^^He film, J. Phys.: Condens. Matter 1, 5263-5274 (1989).

22. A.P. Mackenzie and Y. Maeno, The superconductivity of Sr2Ru04 and the physics of spin-
triplet pairing. Rev. Mod. Phys. 75, 657-712 (2003).

23. G.E. Volovik, Topological invariants for Standard Model; from semi-metal to topological
insulator, JETP Lett. 91, 55-61 (2010); |arXiv;0912.0502|

24. A.A. Abrikosov and S.D. Beneslavskii, Possible existence of substances intermediate be-
tween metals and dielectrics, Sov. Phys. JETP 32, 699 (1971).

25. A.A. Abrikosov, Quantum magnetoresistance, Phys. Rev. B 58, 2788 (1998).

26. A.A. Burkov and L. Balents, Weyl semimetal in a topological insulator multilayer, Phys.
Rev. Lett. 107, 127205 (2011); A.A. Burkov, M.D. Hook, L. Balents, Topological nodal
semimetals, Phys. Rev. B 84, 235126 (201 1).

27. Xiangang Wan, A.M. Turner, A. Vishwanath and S.Y. Savrasov, Topological semimetal and
Fermi-arc surface states in the electronic structure of pyrochlore iridates, Phys. Rev. B 83,
205101 (2011).

28. S. Ryu and Y. Hatsugai, Topological origin of zero-energy edge states in particle-hole sym-
metric systems, Phys. Rev. Lett. 89, 077002 (2002).

29. J.L. Manes, F. Guinea and M.A.H. Vozmediano, Existence and topological stability of Fermi
points in multilayered graphene, Phys. Rev. B 75, 155424 (2007).

30. M. A. H. Vozmediano, M. I. Katsnelson, F. Guinea, Gauge fields in graphene, Physics Re-
ports 496, 109 (2010).

31. A. Cortijo, F. Guinea, M.A.H. Vozmediano Geometrical and topological aspects of graphene
and related materials, arXiv: 11 12.20541

32. A.P. Schnyder and S. Ryu, Topological phases and flat surface bands in superconductors
without inversion symmetry, arXiv; 1011. 1438 , Phys. Rev. B 84, 060504(R) (2011).

33. Andreas P. Schnyder, P. M. R. Brydon, Carsten Timm, Types of topological surface states in
nodal noncentrosymmetric superconductors, arXiv;ll 11.1207

34. A.M. Essin and V. Gurarie, Bulk-boundary correspondence of topological insulators from
their Green's functions, Phys. Rev. B 84, 125132 (2011).

35. M. A. Zubkov, Generalized unparticles, zeros of the Green function, and momentum space
topology of the lattice model with overlap fermions, arXiv: 1202.2524

36. M.V. Berry, Quantal phase factors accompanying adiabatic changes, Proc. R. Soc. Lond. A
392,45-57 (1984).

37. G.E. Volovik, Zeros in the fermion spectrum in superfluid systems as diabolical points,
Pis'ma ZhETF 46, 81-84 (1987); JETP Lett. 46, 98-102 (1987).

38. PG. Grinevich and G.E. Volovik, Topology of gap nodes in superfluid ^^He: K4 homotopy
group for ^He - B disclination, J. Low Temp. Phys. 72, 371-380 (1988).

39. G.E. Volovik, Gapless fermionic excitations on the quantized vortices in superfluids and
superconductors, JETP Lett. 49, 391-395 (1989).

40. M.A. Silaev and G.E. Volovik, Topological superfluid ■'He-B; fermion zero modes on inter-
faces and in the vortex core, J. Low Temp. Phys. 161, 460-473 (2010); I^Xrv:1005.4672l

41. V.A. Khodel and VR. Shaginyan, Superfluidity in system with fermion condensate, JETP
Lett. 51, 553 (1990).

42. G.E. Volovik, A new class of normal Femii liquids, JETP Lett. 53, 222 (1991).

43. Y. Tsutsumi, M. Ichioka, and K. Machida, Majorana surface states of superfluid ■'He A and
B phases in a slab, Phys. Rev. B 83, 094510 (2011).

46

G.E. Volovik

44. G.E. Volovik, Flat band in the core of topological defects: bulk-vortex correspondence in
topological superfluids with Fermi points, JETP Lett. 93, 66 (2011).

45. CD. Froggatt and H.B. Nielsen, Origin of Symmetry, World Scientific, Singapore, 1991.

46. J. von Neumann und E.P. Wigner, Uber das Verhalten von Eigenwerten bei adiabatischen
Prozessen, Phys. Zeit. 30, 467^70 (1929).

47. S.P. Novikov, Magnetic Bloch functions and vector bundles. Typical dispersion laws and
their quantum numbers, Sov. Math., Dokl. 23, 298-303 (1981).

48. H. Georgi, Another odd thing about unparticle physics, Phys. Lett. B650, 275-278 (2007);
arXiv:0704.2457

49. M. Luo and G. Zhu, Some phenomenologies of unparticle physics. Phys. Lett. B 659, 341
(2008).

50. M. Pepea and U.J. Wieseb, Exceptional deconfinement in G(2) gauge theory, Nucl. Phys.
B768, 21-37 (2007).

51. M. Kadastik, K. Kannike and M. Raidal, Dark Matter as the signal of Grand Unification,
Phys. Rev. D 80, 085020 (2009).

52. S. Adler, Axial-vector vertex in spinor electrodynamics, Phys. Rev. 177, 2426-2438 (1969).

53. J.S. Bell and R. Jackiw, A PCAC puzzle: Tt(, — > 77 in the a model, Nuovo Cim. A 60, 47-61
(1969).

54. L. Canetti, M. Drewes and M. Shaposhnikov, Matter and Antimatter in the Universe,
arXiv:1204.4186|

55. T.D.C. Bevan, A.J. Manninen, J.B. Cook, J.R. Hook, H.E. Hall, T. Vachaspati and G.E.
Volovik, Momentum creation by vortices in superfluid ■'He as a model of primordial baryo-
genesis. Nature 386, 689-692 (1997).

56. F.R. Klinkhamer and G.E. Volovik, Emergent CPT violation from the splitting of Fermi
points. Int. J. Mod. Phys. A 20, 2795-2812 (2005); hep-th/0403037

57. L. Tarruell, D. Greif, T. Uehlinger, G. Jotzu and T. Esslinger, Creating, moving and merging
Dirac points with a Fermi gas in a tunable honeycomb lattice. ,arXiv: 111 1 .5020

58. F.R. Klinkhamer, Possible new source of T and CP violation in neutrino oscillations, Phys.
Rev. D 73, 057301 (2006).

59. G.E. Volovik and L.P. Gorkov, Superconductivity classes in the heavy fermion systems,
JETP 61, 843-854 (1985).

60. D.B. Kaplan and Sichun Sun, Spacetime as a topological insulator: Mechanism for the origin
of the fermion generations, Phys. Rev. Lett. 108, 181807 (2012).

61. G.E. Volovik and VA. Konyshev, Properties of the superfluid systems with multiple zeros
in fermion spectrum, JETP Lett. 47, 250-254 (1988).

62. P. Hofava, Spectral dimension of the Universe in quantum gravity at a Lifshitz point, Phys.
Rev. Lett. 102, 161301 (2009).

63. P Hofava, Quantum gravity at a Lifshitz point, Phys. Rev. D 79, 084008 (2009)^

64. P Hofava, Membranes at quantum criticality, JHEP 0903, 020 (2009), arXiv:0 812.4287l

65. Cenke Xu and P. Hofava, Emergent gravity at a Lifshitz point from a Bose liquid on the
lattice, Phys. Rev. D 81, 104033 (2010).

66. E.M. Lifshitz, On the theory of second-order phase transitions I, II, Zh. Eksp. Teor. Fiz. 11,
255, 269 (1941).

67. G.E. Volovik, Reentrant violation of special relativity in the low-energy comer, JETP Lett.
73, 162-165 (2001); hep-ph/0101286

68. P. Dietl, F. Piechon and G. Montambaux, New magnetic field dependence of Landau levels in
a graphenelike structure, Phys. Rev. Lett. 100, 236405 (2008); G. Montambaux, F. Piechon,
J.-N. Fuchs, and M.O. Goerbig, A universal Hamiltonian for motion and merging of Dirac
points in a two-dimensional crystal, Eur. Phys. J. B 72, 509-520 (2009); arXiv:0907.0500]
R. deGail, J.-N. Fuchs, M.O. Goerbig, F. Piechon and G. Montambaux, Manipulation of
Dirac points in graphene-like crystals, Physica B 407, 1948-1952 (2012).

69. YD. Chong, X.G. Wen and M. Soljacic, Effective theory of quadratic degeneracies, Phys.
Rev. B 77, 235125 (2008).

70. S. Banerjee, R. R. Singh, V. Pardo and W. E. Pickett, Tight-binding modeling and low-energy
behavior of the semi-Dirac point, Phys. Rev. Lett. 103, 016402 (2009).

Topology of quantum vacuum

47

71. K. Sun, H. Yao, E. Fradkin and S.A. Kivelson, Topological insulators and nematic phases
from spontaneous symmetry breaking in 2D Fermi systems with a quadratic band crossing,
Phys. Rev. Lett. 103, 046811 (2009).

72. L. Fu, Topological crystalline insulators, Phys. Rev. Lett. 106, 106802 (201 1).

73. X.G. Wen and A. Zee, Gapless fermions and quantum order, Phys. Rev. B 66, 235 110 (2002).

74. B. Beri, Topologically stable gapless phases of time-reversal invariant superconductors,
Phys. Rev. Lett. 103, 016402 (2009).

75. E. McCann and V.I. Falko, Landau-level degeneracy and quantum Hall effect in a graphite
bilayer, PRE 96, 086805 (2006).

76. M. Koshino and T. Ando, Transport in bilayer graphene: Calculations within a self-
consistent Bom approximation, Phys. Rev. B 73, 245403 (2006).

77. F.R. Klinkhamer and G.E.Volovik, Superluminal neutrino and spontaneous breaking of
Lorentz invariance, Pis'ma ZhETF 94, 731-733 (201 1); arXiv: 1 109.6624

78. F.R. Klinkhamer, OPREA's superluminal muon-neutrino velocity and a Fermi-point-
splitting model of Lorentz violation, ar Xiv: 1109.567 II

79. F.R. Klinkhamer, Superluminal neutrino, flavor, and relativitv. larXiv: 1 1 10.2146 1

80. T.T. Heikkila and G.E. Volovik, Fermions with cubic and quartic spectrum, Pis'ma ZhETF

92, 751-756 (2010); JETP Lett. 92, 681-686 (2010); arXiv: 1010.03931

81. F. Guinea, A.H. Castro Neto and N.M.R. Peres, Electronic states and Landau levels in
graphene stacks, Phys. Rev. B 73, 245426 (2006).

82. Kin Fai Mak, Jie Shan and TF. Heinz, Electronic structure of few-layer graphene: Exper-
imental demonstration of strong dependence on stacking sequence. Phys. Rev. Lett. 104,
176404 (2010).

83. A.H. Castro Neto, F. Guinea, N.M.R. Peres, K.S. Novoselov and A.K. Geim, The electronic
properties of graphene. Rev. Mod. Phys. 81, 109-162 (2009).

84. T.T. Heikkila and G.E. Volovik, Dimensional crossover in topological matter: Evolution of
the multiple Dirac point in the layered system to the flat band on the surface, Pis'ma ZhETF

93, 63-68 (2011); JETP Lett. 93, 59-65 (20in; I^Xw:1011.4185l

85. G. Altarelli and F. Feruglio, Discrete flavor symmetries and models of neutrino mixing. Rev.
Mod. Phys. 82, 2701-2729 (2010).

86. T.T. Heikkila, N.B. Kopnin and G.E. Volovik, Flat bands in topological media, Pis'ma
ZhETF 94, 252-258 (2011); JETP Lett. 94, 233-239(2011); arXiv: 1012.0905

87. Jens O. Andersen and Tor Haugset, Magnetization in (2-l-l)-dimensional QED at finite tem-
perature and density, Phys. Rev. D 51, 3073-3080 (1995).

88. M.I. Katsnelson and G.E. Volovik, Quantum electrodynamics with anisotropic scaling:
Heisenberg-Euler action and Schwinger pair production in the bilayer graphene, JETP Lett.
95, 411^15 (2012); arXiv: 1203.1578i

89. M.I. Katsnelson, G.E.Volovik and M.A.Zubkov, Euler-Heisenberg effective action and mag-
netoelectric effect in multilayer graphene, ar Xiv:1206.3973l

90. M.A. Zubkov, Schwinger pair creation in multilayer graphene, Pis'ma ZhETF 95, 540-543
(2012); ai-Xiv: 1204.0138

91. C.L. Kane and E. Mele, Z2 topological order and the quantum spin Hall effect, Phys. Rev.
Lett. 95, 146802 (2005).

92. B.A. Volkov, A. A. Gorbatsevich, Yu.V. Kopaev and VV Tugushev, Macroscopic current
states in crystals, JETP 54, 391-397 (1981).

93. B.A. Volkov and O.A. Pankratov, Two-dimensional massless electrons in an inverted con-
tact, JETP Lett. 42, 178-181 (1985).

94. G.E. Volovik, Fermion zero modes at the boundary of superfluid ^^He-B, Pis'ma ZhETF 90,
440-442 (2009); JETP Lett. 90, 398^01 (2009); arXiv:0907.5389

95. G.E. Volovik, Exotic properties of superfluid ^He, World Scientific, Singapore, 1992.

96. G.E. Volovik, Fractional statistics and analogs of quantum Hall effect in superfluid ■'He
films. In: Quantum Fluids and Solids - 1989 ed. by G.G.Ihas, Y.Takano (AIP Conference
Proceedings , 1989) 194, pp. 136-146.

97. VP Mineev and G.E. Volovik, Planar and linear solitons in superfluid "'He, Phys. Rev. B 18,
3197-3203 (1978).

48

G.E. Volovik

98. M.M. Salomaa, G.E. Volovik, Quantized Vortices in superfluid He, Rev. Mod. Phys. 59,
533-613 (1987).

99. G.E. Volovik, Topological invariant for superfluid ^He-B and quantum phase transitions,
Pis'ma ZhETF 90, 639-643 (2009); JET? Lett. 90, 587-591 (2009); arXiv:0909.3084

100. G.E. Volovik, Topological superfluid ^He-B in magnetic field and Ising variable, JETP Lett.
91, 201-205 (2010); arXiv: 1001.15"T4 ,

101. V. Gurarie and L. Radzihovsky, Resonantly-paired fermionic superfluids, Ann. Phys. 322,
2-119 (2007).

102. R. Jackiw' and C. Rebbi, Solitons with fermion number 1/2, Phys. Rev. D 13, 3398-3409
(1976).

103. Y. Nishida, Is a color superconductor topological? Phys. Rev. D 81, 074004 (2010).

104. T. Ohsaku, BCS and generalized BCS superconductivity in relativistic quantum field theory:
Formulation. Phys. Rev. B 65, 024512 (2001).

105. Y. Nishida, L. Santos and C. Chamon, Topological superconductors as nonrelativistic limits
of Jackiw-Rossi and Jackiw-Rebbi models, arXiv: 1007.2201 II

106. S.M.A. Rombouts, J. Dukelsky and G. Ortiz, Quantum phase diagram of the integrable p, +
ipy fermionic superfluid, ar Xiv: 1008.3 406

107. R. Jackiw and P. Rossi, Zero modes of the vortex-fermion system, Nucl. Phys. B 190, 681-
691 (1981).

108. G.E. Volovik, Localized fermions on quantized vortices in superfluid ^He-B, /. Phys.: Con-
dens. Matters, 357-368 (1991).

109. J.C.Y. Teo and C.L. Kane, Majorana fermions and non-Abelian statistics in three dimen-
sions, Phys. Rev. Lett. 104, 046401 (2010).

1 10. J.C.Y. Teo and C.L. Kane, Topological defects and gapless modes in insulators and super-
conductors, Phys. Rev. B 82, 115120 (2010).

111. Chi-Ken Lu and LF. Herbut, Pairing symmetry and vortex zero-mode for superconducting
Dirac fermions, Phys. Rev. B 82, 144505 (2010)

1 12. C. Caroli, P G. de Gennes and J. Matricon, Phys. Lett. 9, 307 (1964).

113. G.E. Volovik, Vortex motion in fermi superfluids and Callan-Harvey effect, JETP Lett. 57,
244-248 (1993).

114. T. Sh. Misirpashaev and G.E. Volovik, Fermion zero modes in symmetric vortices in super-
fluid ^He, Physica B, 210, 338-346 (1995).

115. T. Mizushima and K. Machida, Vortex structures and zero-energy states in the BCS-to-BEC
evolution of p-wave resonant Fermi gases, Phys. Rev. A 81, 053605 (2010).

116. R.M. Lutchyn, J.D. Sau and S. Das Sarma, Majorana fermions and a topological phase
transition in semiconductor-superconductor heterostructures, Phys. Rev. Lett. 105, 077001
(2010).

117. Y. Oreg, G. Refael and F. von Oppen, Helical liquids and Majorana bound states in quantum
wires, Phys. Rev. Lett. 105, 177002 (2010)

118. C.W.J. Beenakker, Search for Majorana fermions in superconductors, arXiv:1112.1950

1 19. Ken Shiozaki and Satoshi Fujimoto, Green's function method for line defects and gapless
modes in topological insulators: beyond semiclassical appoach. larXiv:1111.1685l

120. Zhong Wang, Xiao-Liang Qi, Shou-Cheng Zhang, General theory of interacting topological
insulators, arXiv: 1004.4229

121. D.B. Kaplan, Method for simulating chiral fermions on the lattice, Phys. Lett. B 288, 342-
347 (1992); arXiv:hep-lat/9206013

122. M.F.L. Golterman, K. Jansen and D.B. Kaplan, Chem-Simons currents and chiral fermions
on the lattice, Phys.Lett. B 301, 219-223 (1993): arXiv: hep-lat/9209003

123. M.A. Zubkov and G.E. Volovik, Topological invariants for the 4D systems with mass gap.
Nuclear Physics B 860, Issue 2, 295-309 (2012); arXiv: 1201.4185

124. N.B. Kopnin and M.M. Salomaa, Mutual friction in superfluid ■'He: Effects of bound states
in the vortex core, Phys. Rev. B 44, 9667-9677 (1991).

125. K. Akama, An attempt at pregeometry - gravity with composite metric. Prog. Theor. Phys.
60, 1900(1978).

Topology of quantum vacuum

49

126. G.E. Volovik, Superfluid ^He-B and gravity, Physica B 162, 222 (1990).

127. C. Wetterich, Gravity from spinors, Phys. Rev. D 70, 105004 (2004); Universality of geom-
etry, Phys. Lett. B 712, 126-131 (2012).

128. D. Diakonov, Towards lattice-regularized quantum gravity, [arXiv: 1 109.00911